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^ ■ Abstract. A LG-WKB and Turning point theory is developed for three term re- 

O , currence formulas associated with monotonic recurrence coefficients. This is used 

^ I to find strong asymptotics for certain classical orthogonal polynomials including 

Wilson polynomials. 
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Ph' I. Introduction 

■ There has been much recent interest in so-called e-difference equations. As shown by 
Deift and and McLaughlin [11] the recurrence formula satisfied by polynomials orthogonal 

^ . on the real, 

■ a„+i?/'„+i(x) + (6^ - x)V'n(a:) +ani/'n-i(a;) = (1.1 

00 ■ 



with a„ > 0, and real can be recast to satisfy such an equation. To see this we assume 
that On, and hn are discretizations of the function a{u) and h{u) and in order to use Taylor 
series we will write u = where in this article t will be restricted to a compact interval 
of the real line. In the examples considered below an auxiliary scaling x = X^y + Xl first 
^ ! suggest by Nevai and Dehesa [33] and developed by Van Assche [43] will be performed 
where and are chosen so that a(t, e) = a(|)/Ae and e) = — A^)/Ae are 

bounded functions of t and e for (t, e) in a compact set say [0, T] x [0, eo]. In most of the 
cases we will discuss Ae = a(j + |) and A^ = 0. This produces the e difference equation, 

a{tn + e, e)ip{tn + e, y, e) + e) - y)ip{tn, y, e) + a(t„, e)?/5(t„ - e, y, e) = 0, (1.2) 

where = ne, and il){tmy-,e) = ?/'n(Aey). Costin and Costin [7] considered higher order 
difference epsilon equations. (For an application of this theory to knots see [13].) In both 
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cases it was assumed that the coefficients in the difference equation were C°° functions of 
t. In the above works the authors developed a LG-WKB theory in the region where the 
roots of the eikonal equation (see e.g. [11]) are distinct. Furthermore they also produced 
solutions in terms of Airy functions in a neighborhood where two roots of the eikonal 
equation collide at a rate (Ai(t) — A2(t)) ~ ^/{t — tp). Here the point of collision tp is 
called a turning point. The neighborhood of tp in which the above solutions are valid was 
shown to be of order \t — tp\ < e^^'^ which is sufficiently large to allow matching of the 
WKB solutions of the "exterior" region to those of the "interior" region around tp. In a 
later paper Geronimo et al [14] developed a discrete analog of the Langer transformation 
used in differential equations. With this transformation the authors gave conditions on the 
recurrence coefficients that allowed solutions, in terms of Airy functions, valid in an 0(1) 
neighborhood of tp. They then went on to use these solutions to find strong asymptotics 
for Hermite polynomials that resembled the well known Plancheral-Rotach asymptotics 
obtained by other means (see Olver [35], Szego [41]). In the process of finding useful 
approximate solutions the authors considered the singular initial value problem. This 
problem arises when the off-diagonal coefficients in the recurrence coefficients, a{tn,e), 
equal or tend to zero with e in a neighborhood of t^ = 0. Such a phenomena occurs 
commonly for instance in the continuum limit of the Toda lattice [11], in the case where 
the recurrence coefficients tend to infinity with n and in many cases in the theory of varying 
recurrence coefficients [25] and [1]. 

The strong asymptotics of polynomials orthogonal with respect to exponential weights 
(i.e., Hermite polynomials, Freud weights, etc.) has received much attention in recent years, 
Baik et al [2], Bleher and Its [3], Deift et al [10], Levin and Lubinsky [28], McLaughlin and 
Kriecherbauer [23], VanLessen [43] (see also Nevai [31, 32], Plancherel and Rotach [34], 
and Sheen [40]). In the above studies the main starting point is in weight function with 
respect to which the polynomials are orthogonal. In these studies heavy use is made of 
the theory of external fields both constrained (Buyarov and Rakhmanov [4] , Dragnev and 
Saff [9], Kuijlaars and Van Assche [26], Lubinsky and Levin [29], Rakhmanov [36, 37]) and 
unconstrained (Gonchar and Rakhmanov [19], Saff and Totik [39 and references therein]). 
Other recent methods used to study these types of polynomials have been developed by Jin 
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and Wong [21], Maejima and Van Assche [27]. Strong asymptotics with error bounds for 
orthogonal polynomials having recurrence coefficients that tend in magnitude to infinity 
have been obtained by Geronimo, Bruno, and VanAssche [14], Geronimo and Smith [16], 
and Geronimo, Smith and Van Assche [17], Geronimo and Van Assche [18], Van Assche 
and Geronimo [42] and Wang and Wong [44] (for a heuristic approach see Dominici [8]). 
In Wang and Wong second order difference equations of the form 

Pn+lix) - (AnX + Bn)Pn{x) + Pn-l{x) = 

are considered where 

6=0 n=0 

For real x approximate solutions in terms of Airy functions having a complete asymptotic 
expansions are found. Use of their technique to study the asymptotics of orthogonal 
polynomials however requires that the solution of the initial value problem be computed 
by other means. Also in this vein are the uncontracted asymptotics obtained by Janus and 
Naboko [20] . It is perhaps worthwhile to emphasize that the theory developed below allows 
us, in a self contained manner, to obtain the strong asymptotics of solutions of initial value 
problems associated with epsilon difference equations along with error bounds. The main 
results of this paper are to obtain strong asymptotics of continuous dual Hahn, Hermite, 
Laguerre, Meixner, Meixner-PoUaczek and Wilson polynomials pn with an error of 0(e) 
where e = 1/n. 

We begin in Section II by reviewing the results of [14] and applying them to monotonic 
recurrence coefficients. Then in Section III we review the theory of external fields and show 
the connection between the discrete analog of the Langer transformation and potentials 
arising from these fields. These results are applied to investigate monotonic recurrence 
coefficients that allow an analytic continuation to a wedge. This allows us to simplify 
formulas for the external field and the Langer Transformation. In Section IV we discuss the 
singular initial value problem which allows us to match solutions of the initial value problem 
with those obtained in terms of Airy functions. The matching problem is considered in 
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Section V where the case of discrete weights is also considered. In Section VI we consider 
some model recurrence coefficients associated with symmetric, antisymmetric, and discrete 
Freud weights. In Section VII special cases of these coefficients are chosen so as to well 
approximate the recurrence coefficients of the classical orthogonal polynomials mentioned 
above. The location of their zeros is also discussed. 

II. Preliminaries 

Tarning Point Theory 
(1) Real case 

We begin by reviewing the results of [14] that will be needed for the rest of the paper. 
The objective is to obtain solutions to (1.2) in terms of Airy functions. We suppose that 
y is a real variable, and in order to exploit the symmetry in equation (1.2) we expand 



a{t + e, e) = a{t + e/2, e) + a{t + e/2, e)'e/2 + 0{e^) and a(i, e) = a{t + e/2, e) - a{t + 
e/2,e)'e/2 + 0(e2). Set 



cosh kit, y, e) 



y - b{t, e) 



for t<t 



cos k{t, y, e) 



2a (t+ f,e) 
y - hit, e) 



for t>t 



(2.1) 



2a (t+ f,e) 



'P-i 





for t < t. 



for t > tp. 



(2.2) 



With 




(2.3) 



we set 



'^i{t, y, e) = g{t, y, e)Ai(e 3p(t, e)), 



(2.4) 
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and 

V'2(i, y, e) = git, y, e)Bi(e-ip(t, y, e)), (2.5) 
where Ai and Bi are Airy functions [35] and satisfy the differential equation 

X{x)" = xx{x). 

From (2.1) it foUows that, 

and equation (2.2) can be written in the compact form, 

2 3 

-p2(t,y,e) = k{u,y,e)du. (2.7) 

Remark. For real variables the above expression should be considered as a short hand 
notation for (2.2). 

Remark. The expansion about t + 1/2 discussed above is motivated by the special 
functions we will consider below and it gives a simpler error bound for the approximation 
problem discussed below. 

The rest of this section is devoted to showing that (2.4) and (2.5) can be used as useful 
approximate solution to (1.2)lJnless otherwise stated the branch of lii{z) which is positive 
for 2; > 1 and analytic in a neighborhood of the positive real axis will be used. In what is 
to follow much use will be made of the analytic properties of the function, 

4y) = y + Vy' - 1, (2.8) 

where the branch of the square root function is chosen so that for large y, z ^ 2y. The 
above function maps the exterior of [—1, 1] to the exterior of the unit circle so that, 

1^1 >l, yeC\[-l,l], 1^1 = 1, ye [-1,1]. (2.9) 

For y real the functions z+ = lim/i_>o z{y + ih) and Z- = lim/i_>o z{y — ih) are continuous 
functions of y. 

We will denote the interior (relative to the complex plane) of a set K as int(i^), 
C+ = {y : Imy > 0}, C_ = {y : Imy < 0} and C± = C± U R. 
We begin with, 
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Lemma 2.1. Let f{z) = \ia^{z + \J z'^ — 1) then f is analytic for z e C\(— oo, —1]. In this 
region f has a sole zero which is simple and located at z = 1. 



Proof. From the mapping properties ot z + \J z"^ — \ we see that In^ {z + \J z^ — is 
analytic for z e C\(— oo,l]. Furthermore Im(ln^ (2; + V-s^ — l)) = for 2; e [—1,1]. 
Consequently by the Schwarz reflection principal In^ {z + V-z^ — l) is analytic for z e 
C\(— 00, —1]. If we set = 2; — 1 and use the representation 



lnfl + t(; + yty(tf; + 2)) = I ^ . 
V /Jo ^x{x + 2) 



we find 



for \w\ < 2. Thus / has a simple zero at 2; = 1. That this is the only zero in C\(— 00, 1] 
follows from the mapping properties of 2; + \/z'^ — 1. □ 
We will suppose that there exists T > 1 and eo sufficiently small so that 

a{t, e),h{t, t) e C([0, T] x [0, e]) (2.11) 

^,%i).C((0.rlx[0.e„l) i = l,...4 (2.12) 

a(0,0) = = 6(0,0) buta(t,e)>0 (t, e) e (0, T] x [0, eo]. (2.13) 

Set 

7±(t, e) = h{t, e) ± 2a{t + e/2, e), 

and when they exist 

Note that from the assumptions on a(t, e) and 6(t, e), 7^(0, 0) = 0. 
We will assume that one of the following cases occur: 

1) For each e G [0, cq] 7'^(t, e), and — 7^(t, e) are increasing for t G [0,T] with nonzero 
derivatives in t for (t, e) G (0, T) x [0, eo]. For every interval [1/1,2/2] G (0, 7o~(T)) there is 
an ei such that t+ G C([j/i, ^2] x [0, ei]). Likewise for every interval [2/1,^/2] £ (7(7 (2^)) 0) 
there is an e\ such that t~ G C([yi,y2] x [0,ei]). 



2) For each e e [0, eo] , 7"*" (t, e) is increasing and 7~ (t, e) = for t e [0, T] . The derivative 
of 7"'"(t, e) with respect to t is nonzero for (t, e) e (0, T) x [0, eo]. For every interval 
[1/1,1/2] e (0,7^(T)) there is an ei such that t+ G C([yi,y2] x [0, ei]). 

3) 7+(t,e) and 7"(t,e) are increasing for t G [0, T] with nonzero derivatives in t for 
(t, e) e (0, T) X [0,eo]. For every interval [2/1,2/2] € (0, 7o^(2^)) there is an ei such that 

e C([yi,y2] X [0,ei]). Likewise for every interval [2/1,^/2] £ (0,7o'(^)) there is an ei 
such that t+ G C{[yi,y2] x [0, ei]). 

We restrict our attention to the above cases since these are the ones that are relevant 
for the examples we wish to consider. Other cases for instance ^^{t, e) and 7~(t, e) both 
decreasing or j~^{t, e) = and 7~(t, e) decreasing are also important and can be treated by 
interchanging the roles of 7+ and 7". Set 7^(t) = 7"''(t, e), 77 (t) = 7~(t, e), z'^{t,y,e) = 
and z-(t,y,e) = 7rTr#7^- Note that z+ - 1 = /,T?;%? , and z' - 1 = 

2a(t+e/2,e) ' V / 2a(t+e/2,e) 2a(t+e/2,e) 

2a{t+e/2,e) " 

Lemma 2.2. Suppose that a{t, e) and e) satisfy conditions (2.11)-(2.13) and case 1, 2, 
or 3 holds. Then for every [tin, tfi] x [yi, j/2] C (0, T) x (0, 70" (T)) there is an ei suci that 

& ^'^+lye)-t ^ '^([tin^tfi] X [1/1,1/2] X [0, 61] ) for z = 0, . . . , 3. In case 1 and y < for every 
[t^n,tf^] X [2/1,2/2] C (0, T) X (7o"(7^), 0) there is an ei such that e C([tin, tfi] x 

[2/1,2/2] X [0, ei]) for z = 0, . . . , 3. Finally for case 3 tiiere exits an ei suci tiat for every 
interval[tin,tfi]x[yi,y2] C (0, T) x (0, 7o-(T)), 9L5Z^h^^ g C([tin, t/i] x [yi, ^2] x [0, ei]) 
for z = 0, . . . , 3. 

Proof. Fix [2/1,2/2] £ (0,7o~(^)) by enlarging [/:in,^/i] if need be and using the continuity 
of 7(J" we can assume that [2/1, 2/2] C 7(J"((tm, ^/i))- Choose ei so that t^{y, e) G C( [2/1, 2/2] x 
[0, ei]) and [2/1,2/2] C 7e^((^m, ^/i)) for aU e G [0, ei]. For {y, t, e) G [tin, tfi] x [2/1,2/2] x [0, ei] 
and z = 0, . . . , 3 set, 

^/(t 2/ e) - I - '^'^^'^'^ 

Q^iJKhy,^) S ^-^+fye^ 

To show that / G C {[tin, tfi] x [2/1,2/2] x [0,ei]) it is sufficient to demonstrate that the 
above function is continuous at t = t^{y,e). Other values of {t,y,e) follow from the 



the smoothness conditions imposed on 7+ and a{t,e). Fix (t = tp{y,e),y,e) and 
suppose (t, y, e) is a point close by. For (y, e) close to (y, e), (y, e) is close to t^{y, e). If 
t = tp{y, e) then use the second part of the definition of / otherwise expand 'j~^{t, e) in a 
Taylor polynomial with remainder about tp{y, e). These formulas give the contiuity of /. 
An argument similar to that above shows that ^/(t, y, e) e C{[tin-tfi] x [yi,j/2] x [0,ei]), 
z = 0, . . . , 3. An analogous discussion for case 1 with y < and the second part of case 
three gives the result. □ 
With the above we can now discuss the smoothness properties of k^. 

Lemma 2.3. Suppose that a{t, e) and b{t, e) satisfy conditions (2.11)-(2.13) and case 1 or 
2 holds. Then for every [tin,tfi] x [2/1,^/2] C (0,r) x (0,7q"(T)) there is an ei such that for 

k given by equation (2.6), 

& tUy.T-t ^ ^(t^^'^' ^^'^l [^1' x [0, ei]), i = 0, . . . , 3 with strictly positive 

in [tin, t ft] X [yi,y2\ x [0,ei], and 
.. > 0, (l^^f^l^l > 0), for {y,t,e) e x [y„y,] x [0, 61] and 

t<t+iy,e) it>t+iy,e)). 
For case 1 and y < set 



k\t, y, e) ^In [ ^ ^' x + 1/ ( ^ 7'^' ^' J . 1 (2-14) 
^'^'^ ' 2a(t + e/2,e) V V2a(t + e/2,e) ' ' ^ ^ 

Then for every [tin,tfi] x [j/i, 7/2] C (0, T) x (7^(T), 0) tiere is an ei so tiat ij, and ii) hold 
with tp replaced by t~ . 

Proof. Fix [2/1,2/2] e (0,7(|(T)) and choose ei so that t^{y,e) e C([2/i,2/2] x [0, ei]). For 
(t, 2/, e) e [tin, tfi] X [yi, ^2] X [0, ei] and z = 0, . . . , 3 set, 



dPt+{y,e)-t ^T^pKV,^) 



Lemmas 2.1, 2.2, the conditions on 7+ and its derivative with respect to t show that 
h is positive and -^h e C{[tin,tfi] x [2/1,2/2] x [0, ei]), i = 0, . . . , 3. Part ii) of the Lemma 
follows from part i) and equation (2.1). For y < we see from equation (2.14) that 



k'^{t~{y,e),y,e) = thus the argument above can be taken over to this case so that the 
second part of the Lemma follows. □ 

For case 3, 7" becomes an obstacle to the smoothness of k^. Thus the intervals to be 
considered will need to be restricted. 

Lemma 2.4. Suppose that a(t, e) and b{t,€) satisfy conditions (2.11)-(2.13) and case 3 
holds. Then for each [tin.tfi] x [1/1,2/2] G (0,7") x (7(^(T), 7(j~(T)) there is an ei such that 
i), a) of Lemma 2.3 hold. If yi < JgiT) then for each [tin,tfi] x [2/1,2/2] with [2/1,2/2] ^ 
(0, 7q"(T)) and [tin,tfi\ C (0, (7(7)~^(yi)) there is an ei such that i), ii) of Lemma 2.3 hold. 
If 

(t. e) = ( ^ , IfJM. ~ 



k\t, 2/, e) ^ In" ; ' \ [ ; / ^ -1 (2-15) 

^'^'^ ' 2a(t + e/2,e) \l \2a{t + e/2, e) J ' ^ ^ 

then for [Unjfi] x [2/1,2/2] witii [2/1,2/2] C (0,7o"(T)) and [tin,t/i] C ((7(j")"^(2/2), T) there is 
an ei sucii tiiat t~ (2/, e) = {%)~^{y), {y, e) G [2/1, 2/2] x [0, ei] exists and i), ii) of Lemma 2.3 
hold with t+ replaced by t~ and in i) -p— - replaced by For [2/1,2/2] C (0,7^(T)) 

tien it is possible totfi < tin- 

Proof. Supppose that [2/1, 2/2] ^ (0, 7o"(T)) and ei is such that that (2/, e) = (7^)~"'^(2/) ^ 
C{[yi,y2] X [0,ei]). Likewise for any interval [2/1,2/2] G (0, 7o"(^))7 there is an ei such that 
t-{y,e) = (7e")"Hy) e C([2/i,2/2] x [0,ei]). FYom a(t,e) > for (t, e) e (0,T] x [0, eo] 
we see that 7+(t,0) > 7-(t,0) for t e (0,T). Thus for [tin,tfi] x [2/1,2/2] e (0,T) x 
(7q"(T), 7(]~(T)) by (2.11) we can choose ei sufficiently small so that y > %{t) for all 
{t,y,e) G [tin7^/i] X [2/1,2/2] X [0, ei]. That i) and ii) hold now follows from an argument 
similar to that given in Lemma 2.3. Suppose now 2/1 < 1q (T) and [tmj^/i] x [2/1,2/2] is 
such that [2/1,2/2] e (0,7(1(2^)) ^■nd ([tirn^/i] C (0, (7(7)~^(yi))- As above ei may be chosen 
sufficiently small so that y > ^~{t) for all (t, y, e) e [tin, tfi] x [2/1,2/2] x [0, ei] so that i), and 
ii) will be satisfied. Suppose now that [Un, t/i] x [2/1,2/2] is such that [2/1,2/2] C (0, 7~(T, 0)) 
and [iim ifi] C ((7o~)~"'^(2/2), T). Then the uniform continuity of a{t, e) and 6(/:, e) given by 
(2.11) shows that by choosing ei sufficiently small we can insure that y < 7^(t) for all 
{t,y,e) e [tin, ifi] X [2/1,2/2] X [0,ei]. Thus smoothness part of condition i) for k"^ given by 
(2.15) follow from Lemmas 2.1 and 2.2 and an argument analogous to the one used above. 
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Since A;^ > for y < 7^ (t) and has a simple zero at j/ = 7^ (t) the properties of tp (y, e) and 
^77 imply that ^ > Q. Property ii) follows by the choice of [im, ifi] x [yi, 1/2] x [0, ei] 

and the fact that sinh^ /c = — 1. The fact that 70" > 7,^ for t > allows us to 

choose tfi <tin- □ 
For cases 1 and 3 we will have need of 



and 




(2.16) 



(2.17) 



With the above results we will be able to obtain approximate solutions to the difference 
equation. 

Remark. In the proof of the Theorem below we will use Theorem 3.3 of [14] where it is 
assumed that k'^ is monotonic. While this is true in our case examination of the proof of 
Theorem 3.3 shows that only conditions i)-iii) or ia)-iiia) in [14]. Note that in ii) and iiia) 
of [14] there is an error, the running index should be i = 0, 1, 2, 3. 

Theorem 2.5. Suppose that a{t, e) and b{t, e) satisfy conditions (2.11)-(2.13) and case 1, 
or 2 holds. Let pi be given by equation (2.2) and il^iit, y, e) i = 1, 2 be given by equations 
(2.4) and (2.5) respectively. Then for each [tm,^/i] x [^172/2] C (0, T) x (0, 7(^(T)) there is 
an ei such that for all {t, y, e) e [tin, tfi] x [Vij 1/2] x (0, ei], i/ji satisGes, 

a{t + e, e)'0i(t + e, y, e) + a(t, e)'0i(t - e, y, e) 

(2.18) 

- 2a{t + e/2, e) cosh k{t, y, e)ipi{t, y, e) = f3i{t, y,e) i = 1,2 
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where ipi and Pi e C{[tin,tfi] x [^1,^2] x (0,ei]). For {t,y) e [tin,tfi] x [yi,y2], Pi i = 1,2 
satisfy the inequahties, 

\Pi{t,y,e)\<c{y)e^ sup [|Ai(e-ipi(w, y, e))| + |Ai'(e-ipi(«, y, e))|j , (2.19) 

ue{t-e,t+e) L 

and 

\p2it,y,e)\ < c{y)^ sup [|Bi(e-ipi(«, y, e))| + |Bi'(e-ipiKy, e))|l , (2.20) 

with c{y) uniformly bounded on [^1,2/2]- For case 1 with [tinitfi] x [^1,2/2] C (0,T) x 
(7^(T),0) iet p2 be given by (2.16), k'^ by (2.14) and 'ipi{t,y,e) i = 1,2 by equations 
(2.4) and (2.5) respectiveJy with the appropriate substitutions. Then there is an ei such 
that for all {t,y,e) E [tin,tfi] x [1/1,1/2] x (0,ei], ipi and Pi,i = 1,2 satisfy (2.18) with 
il^i, Pi G C{\tin-,tfi] X [1/1,2/2] X (0, ei]) and Pi satisfies the above bound with pi replaced by 
P2- Finally suppose case 3 holds. Then for each [tin,tfi] x [^1,^2] G (0, T) x (7^^ (T), ^q{T)) 
or if yi < 7^(T) then for each [tin,tfi] x [^1,^2] is such that [^1,^2] e (0,7^(T)) and 
([^in, */i] C (0, {^o)~^iyi)) there exists an ei such that for all {t, y, e) e [tin, tfi] x [2/1, t/2] x 
(0, ei], Vi satisfies, (2.18) witii and /3, G C([t,„,t/,] x [2/1,2/2] x (0,ei]). For {t,y) G 
[^in, V»] ^ [2/i;2/2], Pi i = 1,2 Satisfy the inequalities (2.19) or (2.20) respectively with c{y) 
uniformly bounded in [2/1,2/2]- If is given by equation (2.15) and p2 by equation (2.17), 
then for [iin,ifi\x[yi,y2] is such that [^1,^2] C (0, 7-(T, 0)) and [fj^ J/i] C {{^^)-^{y2),T) 
there is an ei such that for all {t, y, e) G [iin,ifi] x [2/1, 2/2] x (0, ei], ipi satisfies, (2.18) with 
pi replaced by p2- The functions if^i and Pi G C {[tin, tfi] x [2/1,2/2] x (0, ei]). Furthermore 
for {t, y) G [tin, tfi] X [2/1,2/2]? A ^ = 1, 2 satisfy the inequalities (2.19) or (2.20) respectiveJj 
with c{y) uniformly bounded in [yi , 2/2] • 

Proof. Suppose 2/ G (0,7+(r,0)) and set p{t,y,e) = ( ff^'ffj for (t,,y,e) G 
[^m,*/i] X [2/1,2/2] X [0,ei] where [tin, tfi] x [2/1,2/2] x [0, ei] are the sets given in the first 
part of Lemma 2.3. Lemma 2.3 shows that p is positive and ■^p{t,y,e) G C([tin,tfi] x 
[2/1, 2/2] X [0, ei]), i = 0, . . . , 3. For (t, , 2/, e) G [t^n, t/i] x [2/1, 2/2] x [0, ei] set 

(t+(i;,e)-t)3/2 ^""''Vp (y, e) - uY''^p{u,y, e)du for t < ^^(2/, e) 
9(i,y,e)=^ ' ' |p(t+(2/,e),2/,e) iort = tp{y,e) (2.21) 

f4iy,e)(^ - *p(y^ e))i/2p(w, 2/, e)d« for t > tp{y, e) 
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Again the issue is to show that q is continuous at {t = t^{y, e), y, e) for fix (j/, e). Lemma 2.3 
(for the definition of p{t'^{y,e),y,e)) and the mean value theorem for integrals (see 
Lemma 3.1 of Olver [35, p. 399]) shows that q{t,y,e) G C{[tin,tfi] x [1/1,1/2] x [0, ei]) 
and is positive. Next set. 



d 



(t+(2/,e)-t)5/2 Jt 



!^^^'''\t+{y,e) - uf/'i^p{u,y,e)du for t < t,{y,e) 



■ft+{y,e)(.'^ ~ ^P^y^ ^)f''^^Pi'^^ ^)du for t > tp{y, e). 



above definition is self consistent and §iq{t,y,e) G C{\tin,tfi\ x [2/1,1/2] x [0,ei]) which is 
positive. Continuing on in this way gives that 4jjq{t,y,e) G C{[tin,tfi] x [2/1,2/2] x [0, ei]) 



Integration by parts, Lemma 2.3, and the mean value theorem for integrals gives that the 

d 

for i = 0, . . . , 3. Since = (^q(t, y, e))^/"^ we obtain that is positive and 

' ' t+(y,e)-t ^2"iV t+{y,e)-t ^ 

'BF t+{yef-t ^ ^([^»n,^/i] X [2/1,2/2] X [0,ei]) for i = 0, ...,3. An analogous argument for 
the sets in the second part of Lemma 2.3 shows that for case 1 and y G (7~(T, 0),0), 
is positive and fr^^g^ g C{[tin, tfi] x [2/1, 2/2] x [0, ei]) for z = 0, . . . , 3. Case 3 
follows from Lemma 2.4 and arguments similar to those given above. For y > that there 
are solutions well approximated by equations (2.4) and (2.5) follows from Theorem 3.3 in 
[14]. Case 1 with y G (7~(T, 0),0) also follow from Theorem 3.3 in [14] after performing 
the substitutions indicated. The proof for case 3 holds in a similar manner after the 
substitutions indicated. □ 

Remark. Because of the connection to the Airy functions ipi, i = 1,2 satisfy the differ- 
ential equation [35, p. 396], 



where is given by equation (2.6) and 

a2(t+ f)sinh^fc(t)^ ^ 



w{t) 

and 



k{ty 
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Set, Uj{x) = g{x)uj{x), j = 1,2 where 

1/2 



Uj{x) 



—X 



'V3 Vs 



(2.22) 



where i^^/g and Hf'J^ are Hankel functions of the 1st and 2nd kind respectively. An 
important property of these functions is that they do not vanish for real x [45] . Consider 
the difference equation, 

ai((n + l)e, e)/(n + 1) + ai(ne, e)/(n - 1) - (y - 6i(ne, e))/(n) = 0, (2.23) 

where for < e < eo, 

sup |a(t,e) - ai(t,e)| = O(e^) = sup |6(t, e) - 6i(t, e)|, (2.24) 

tE[tin,t fi] tE[tin,tfi\ 

for every interval [tini tfi] C (0, T), then the techniques leading to Theorem 3.8 in [14] give 
the following. 

Theorem 2.6. Suppose that the hypotheses of Theorem 2.5 hold, that ai{t,e) and 
bi{t,€) G C([0,T] X [0, eo]), satisfy (2.24) and that ai(t,e) is strictly positive on every 
compact subset of (0, T] x [0, eo]. For cases 1 and 2 with y > or the first part of case 3 let 
[tinitfi] X [yi,y2] X [0,ei] given from Theorem 2.5. Then there exist solutions fi, i = 1,2 
of (2.23) such that for each [y, e) e [^1,^2] x [0, ei] and all ne e [tim tfi], 



fi{n) = 'ipi{n) + ri{n). 



where 



ri{n) 




Ui(n) 





Ui{n) 



< ce. 



(2.25). 



The constant c is uniform on [tin, tfi] x [1/1,2/2] x [0, ei] and G C{[yi,y2] x [0, ei]). 

For case 1 and y < or the second part of case 3 with [tin, tfi] x [2/1,2/2] x [0, ei] given 
in Theorem 2.5 there exist solutions fi, i = 1,2 of (2.23) such that for each (y, e) G 
[yi,y2] X [0,ei] and aline G [tin, tfi], 



ri[n) 




Ui{n) 





Ui{n) 



< ce. 



(2.26) 



where the constant c is uniform on [tin, tfi] x [2/1, 2/2] x [0, ei] and G C([2/i, 2/2] x [0, ei]). 
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Proof. The result for cases 1-3 and y > foUows from Theorem 3.8 in [14]. For case 1 
and J/ < or the second part of case 3 apply Theorem 3.8 of [14] to (— l)'^/i(n). □ 

III. External Fields 

In order to facilitate the matching problem considered below and to consider solutions 
of (1.2) with y complex we will find it convenient to use the theory of potentials with 
external fields. 

Let Q{x) be a continuous function on a closed interval E of the real line. If S is 
unbounded we asume that lim|2.|^oo = +oo, x E Tj. The above assumptions are 
not the most general possible (see [39]) but they are sufficient for the problem we are 
considering. Let M(E) be the collection of positive Borel probability measures on E with 
compact support and define 

Iq[u) = - j j log \z - t\du{z)du{x) + 2 J Q{t)du{t) 

u e M(S), and E = mf{u e S : Iq{u), exists} then there is a unique equilibrium measure 
V such that E = Iq{v) (see Saff-Totik [39], Rakhmanov [36], Gonchar-Rakhmanov [19]). 
Set V^z) = / log j^\dv{t) and Fq = E - JQ dv, then Viz) + Q(z) > Fq on E and 
V'"{z) + Q{z) = Fq on supp((5). The function Q is called the external field. 

In order to consider discrete measures we need to generalize the above problem to 
a constrained equlibrium problem. Let S be a closed bounded intervel and a a positive 
measure supported on this interval with continuous potential and cr(S) > 1. Let M"'(T,) 
be the set of Borel probability measures such that if e M'^(S) then fi < a. It was proved 
by Dragnev and Saff [9, Theorem 2.1] and Rakhmanov [37, Theorem 3] (see also [24]), that 
for Q as above and constraint cr there is a unique Borel probability measure Uq e M'^(S) 
such that for a constant Fq we have 

t/-g+Q<F^on supp(z.5), 
U^-+Q = Iq on supp{iy^) n supp((7 - v^), 
U^^^+Q>Fq on snpp{a-iy^). 

These results have been applied to difference equations by Deift and McLaughlin [11], 
Kuijlaars and VanAssche [25] and Kuijlaars and Rahkmanov [24] and we have the useful 
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Theorem 3.1. [24, Theorem 9.2] Suppose that a{t) > and b{t) are continuous for 

t e [0,T], I \na{u)\du < oo, 7''~(t) = b{t) + 2d{t) has at most one extremum, which if it 
exists is a maximum and 7~(t) = b{t) — 2d{t) has at most one extremum, which if it exists 
is a minimum. Set < c <T 



Q{x) 



f 

Jo 



log 



X — b{u) 



2d{u) ^ y 2d{u) 



X — b{u) 



where t-{x) = inf{0 < t : x e [y (^), 7"'"(^)]} sind 



cr = / UJudu, 

Jo 



with 



dujt{x) 



1 



1 



= a;.. 



dx 7r^(^+(t)-x)(x-7-(t)) 
Then the constrained equihbrium problem with external field Qc = ^ ^nd constraint 



- has the solution 

c 



c Jo 



(3.1) 



As noted by [29] the condition J \ \nd{t)\dt < oo is needed to insure that Q{x) exists. 
We will assume that 7"'"(t) 7^ 7~(t), < t <T which implies that d{t) > 0, < t < T. 
Kuijlaars and Van Assche [25, Remark 1.5] have noted that for < s < c the support of 
the equilibrium measure is in the interval 

inf 7"(s), sup 7+(s) 

0<s<c 0<s<c 

Furthermore from the definition of the constraint we see [24] that it is not active if 7" 
is decreasing on [0, c] and 7+ is increasing over the same region. However if this is not 
the case then the constraint will be active on the interval [7~in' 7~(c)] interval 
[7^(c),7max] where 7";^ = mint^[o,c] and 7+^^ = maxte[o,c] 7+(^) • 

For X e R\supp{wt), set 

(3.2) 



wi — \wi\ 
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Remark. Note that equation (3.1) is to be interpreted as 

1 r 

u,{E) = - / u;u{E)du, (3.3) 
c Jo 

where E is any Borel subset of supp(^'c) and since supp(c<;t) C supp(^'c) this can be rewritten 
as 

I fix)diyc = - I I f{x)dujudu, (3.4) 

where / e C[supp(t'c)] the set continuous function on supp(z^c). 
If 



g\z) =ln 



z - b{t) 



2a{t) W 2a(t) 

then is the complex Green's function associated with the interval [b{t)—2a{t), b{t)+2a{t)], 
< t < p. From the mapping properties of z{y) (see equation (2.8)) we see that is 
analytic for z e C \ {—oo, b{t) + 2a{t)] and g\. = g^\c± is continuous on C±. Furthermore 
for < t < /5 



z - b(t) 



- 1 



(3.5) 



27Ti 



gl-g'_ = < 2m /J ' dwu x G supp(cc;t) 



We also find that, 



gX + gl=< 0, 



X < 7-(t) 

X G supp(o 
0, x> 7+(t). 

2Reg\ x<^-{t) 



X e supp(a;t) 
2Re^*, a;>7+(t). 



(3.6) 



(3.7) 



Remark. It was noticed by [29] that for x > ^~^{t) 

^ f . ds 
Reg' 



7+(t) V(s-7+(t))(^-7-(t)) 



w'fds, 



and in case 1 for a; < 7 (t) 



Reg' 



J X 



ds 



7+(t) 
7-(t) 



TT 



w[ds. 



(3.8) 



(3.9) 



These formulas are easily obtained by differentiating Re g' with respect to x then integrat- 
ing and using equation (3.2). 
Set 

V%z) = - I \n{z-x)du^{x), (3.10) 

for z G C\ (-oo,supo<s<c7+(s)] 
We have the simple, 
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Lemma 3.2. With the hypotheses of Theorem 3.1, V^{z) is analytic for z e C\ (—00, 7+g 
and 



.rcr ^ 1 f\ Z-b{u) Z-h{u) \ 

Viz) = - / In — + A — - 1 
^ ^ cJo y 2a{u) \^ y 2a{u) J 



2 "\ 

du + lc. (3.11) 



Furthermore = V^\c± have continuous extensions to C± which satisfy 



'min 



27ri, X < In 

0, X> 7+a 



Here lc= \ Jq In a (t) (it. 

Proof. The analyticity properties of follow from its definition. For ^ e C \ (—00, b + 
2a\,a > 0,b real, the useful integral representation holds [39], 




ll+lna=-/ lniz-x)^== (3.13 
Jb-2a \/4:a^ -ix-b)^ 



The integrals from < u < c of g± exist by the hypothesis of Theorem 3.1 and give 
rise to continuous functions. Equation (3.11) can now be obtained from (3.4) with / = 
ln(2; — x) and (3.13). The second part of the Theorem is a consequence of the continuity 
and integrability of g^., (3.6), and (3.1). □ 
We will now restrict the coefficients to those cases that arise in the examples to be 
considered. For t G (0, T) assume one of the following statements hold: 

1) 7~(t) is a strictly decreasing function of t and ^~^{t) is a strictly increasing function 
of t. 

2) 7~(t) is a constant and 7"'"(t) is a strictly increasing function of t. 

3) 7~(t) and ^~^{t) are strictly increasing functions of t. 

In all cases we will assume that 7"''(0) = 7~(0). With these conditions Theorem 3.1 
insures that for cases 1-3 and for each 7"'"(T) > x > 7"'~(0) there is an external field. Also 
in case 1, for each 7~(0) > x > 7~(T) there is an associated external field. It is not 
difficult to see from the definition of Q that lim.j;_,oo{Q{x) — Inx) = 00 for the above cases 
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and \imj;^_^{Q{x) — In |a;|) = oo for case 1. The support of the equihbrium measures for 
the cases above are given in Figure 1. 

Case 2 



Case 1 ' y+{t) 




Y+(t) 




Y~(t) 



y-i (t) 



Case 3 



constraint 
active ^" 




Y+(t) 

constraint not active 

Y"(t) 



c Xx 



Figure 1: Support of the equilibrium measures for cases 1-3 
It is not difficult to compute the densities for these measures [4], [29]. 

Lemma 3.3. The density of the equihbrium measure is given by 

c Jf(^) 



where for case one, 



while for case two 



'fix) 

^-{c)<x<0 
(7+)-i(a;), 0<a;<7+(c) 

f{x) = {^+)-\x), 0<x<7+(c). 



For case 3 



where 



du 1 /"-^^^^ 

— ^ = - / Lo'tdt, 0<x< 7+(c) 

dx C 



fix) 



(7~)~^(a;), 0<a;<7"(c) 



c, 



7 (c) < X < 7"'"(c) 



(3.14) 



(3.15) 



(3.16) 
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Proof. The proof follows by substituting the above formulas into the left hand side of 
(3.4) then interchanging the order of integration (see Figure 1) we obtain 



/ g{x)dvc^ / - / g{t)—dx 



where g e C[supp vd- The result now follows from formula (3.4) and the definition of Vc- 

Let ^ be the extension of |^| to x G i?\supp z/^ and ^ be the extension of |^^^^^^^| 
to a; e i?\supp (o" — z^c)- We now have the following, 

Lemma 3.4. Let be as above then, for case 1 

V^{x) + Vl{x) - 2lc - 2Qc{x) = 1 0, 7~(c) < a; < 7+(c) (3.17) 

i-27r c^^>c, 7+(c)<a;<7+(T) 

while for case 2 

{ 0, 7~(c) <x< 7+(c) 

V^{x) + Vl{x)-2l,-2Q,{x)^{ (3.18) 

1^ -27r dvc, 7+ {c)<x< 7+ (T) 

and for case 3 

[ +27r/J"^'=^cZi>e, 0<a;<7-(c) 

V^ix) + Vl{x)-2k-2Q,{x)=l 0, 7-(c)<a;<7+(c) (3-19) 

i-27r dz>e, 7+(c) < a; < 7+(T). 

Proof. The 2nd line of cases 1 and 3 and the first line of case 2 follow from integrating 
the 2nd line of equation (3.7). For the bottom line of the above cases integrate the last 
line of equation (3.7) to obtain 



The right-hand side of the above equation exists from the hypothesis of Theorem 3.1 and 
the continuity of Rg*. Extend the integral to {j~^)~^{x) and use the definition of Q to 
obtain 

2 /■(7+)-'(a;) 
V^ + V1= 2Qc{x) + 24 - - / Rg^dt. 

^ J c 
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Since Re^* = j^^^~^Wt{y)dy (see equation (3.8)) the last line of the above cases follows. 
The first line in case 1 may be derived in an analogous manner by integrating the 1st line 
in (3.7). To obtain the 1st line in case 3 integrate the first line in (3.7) to obtain (see 
Figure 3) 



C 



+ + Reg^dt 

Jiy+)-^x) J{l-)-^{x) 



The first integral is equal to Q while the second is equal to zero (since Re ^r* = 0). The result 
now follows from (3.9) and the definitions of the equilibrium and constraining measures. □ 
We now consider the continuity in e of the above functions. 

Lemma 3.5. Suppose that a{t,e) and b{t,e) satisfy conditions (2.11)-(2.13) and there is 
an integrable m(t) such that 

\lii{a{t,e))\ <m{t), 0<t<T 

for e G [0,eo]. Then for cases 1-3, tK[(y,e) G C(C+ x [0,T] x [0,eo]) while tVl{y,e) G 
C(C_ X [0, T] X [0, cq]). Furthermore for every interval [yi, ^2] C (0, 7(^(T)) there exists an 
ei such that Q G C([yi,j/2] x [0,ei]). For case 1, and y G (7(7(T),0) the same conclusions 
hold with [2/1,^2] C (7^(T),0). In case 1 limj/_>o+ Q{y,0) = = limy_>o_ Q{y,0) while 
in cases 2 and 3 limj/_>o_,_ Q{y, 0) = 0. 



Proof. For fixed e, 7^(t), and a(t, e) satisfy the hypotheses of Theorem 3.1. li y > then 
for cases 1-3 we find from the definition of t- that t_ = tp{y, e). Consider now, 



Q{y,e)= [ 

Jo 



In 



y - h{u,t 



2a{u + f , e 




y - b{u, e) 
2a{u + |, e) 



du. 



For [7/1,^2] C (0,7^(T)) let ei be such that t+(y,e) G C{[yi,y2] x [0,ei]). The inte- 
grability of m as well the fact that the integrand is uniformly continuous for (t, y, e) G 
[^1,^2] X [^1,1/2] X [0, ei] for any [^1,^2] C (0, T) gives the continuity property for Q. From 
equation (3.11) with h{t) = b{t, e) and d{t) = a{t + |, e) for z G C\R we find. 



tV\z) 



g'^{z, e)du + 



In a [u+ 2 ' ^ ' 



(3.20) 
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Equation (3.5) and (2.11) show that the function g*{z, e) is uniformly continuous on com- 
pact subsets of (C\R) x (0, T] x [0, eo], while g±{z, e) are uniformly continuous on compact 
subsets of <C± X (0, T] x [0, eo] respectively. The continuity properties of now follow from 
an argument similar to that given for Q. To show for case 1 that limy_>o_,_ Q{y,0) — 
note that for cases 1-3, (2.13) implies that limy^o+ *p (j/^O) = 0. The integrability of 
now gives the result. Case 1 and y < follows in a similar manner. □ 

Remark. With the assumptions on a{t, e) and b{t, e) above, Lemmas 3.2 and 3.4 are 
easily extended to the case when e > and we will use the the same numbering. 

Prom equations (2.6) and (3.11) we have how p is related to the equilibrium measure. 

Theorem 3.6. Suppose (2.11)-(2.13) hold then for cases 1 and 2 and every rectangle 
[tin,tfi\ X [yi,y2] C (0,T) x (0,7+(T)) or for case 3 with [tin,tfi] x [^1,^2] C (0, T) x 

(7(7 (T), 7+ (T)) or with [t.n.tfil X [^1,^2] such that [^1,^2] C [0,7+(T)] and [^1,^2] C 
(0: (7o')~H2/i)) there exists an ei such that for aU (t,y,e) e [tin,tfi] x [1/1,1/2] x [0, ei], 



•5 J^+{t,e) 

For case 1 and every interval [tin,tfi] x [2/1,1/2] G i^jT) x {^q{T), 0) there exists an ei 
such that for aU (t,y,e) e [tin,tfi] x [2/1,2/2] x [0,ei], 



Finally for case 3 and every interval [ti, x [2/1,2/2] such that [2/1,2/2] € (0, 7o (T")) and 




(3.21) 




(3.22) 



[iin,ifi\ C ((70") ^{y2),T) there exists an ei such that for aU {t,y,e) e [iin,ifi] x [2/1,2/2] x 
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Proof. For cases 1-3, y e [2/1,2/2] and t > t^{y,e) the result follows by using (3.8) 
to rewrite cosh"^ ( 2a{u+e/2e) ) (2-^) ^ integral then interchanging the order of 
integration (see equation (3.8)) using Fubini's theorem. The bottom line of equation (3.22) 
can be obtained in a similar manner (see (3.9)). For the top line in (3.21) apply the above 
procedure to the bottom line of (2.2). The top lines of equation (3.22) and (3.23) follow 
in a similar manner using the bottom lines of (2.16) and (2.17). For the top line of 
equation (3.23) y < ^~(t,e), so recasting cosh~^ as an integral in the top line of (2.17) 
and the definition of i>t gives the result. □ 
(2) Complex case 

We now consider complex extensions of the above results. Motivated by the examples 
discussed earlier we will assume that 

a{t,e) = aq{t,e), b{t,e) = bq(t + ^,ej , a > 0, 6 > 0, (3.24) 

with 

Ut) = q{t,e) = -^^ (3.25) 

Set qe{t) = ^e(t+ 1) then the inverse functions of 7+ and if 6 7^ 2a 77 which we denote 
respectively as t+(y,e) and t-{y,e) are given by t+{y,e) = qe~^{h+2^) and t-{y,e) = 
(le~^ {-h^2^) ■ With Qs = {z ■ |arg(2;)| < 5} we will assume that, 
ic) q{t) is nonegative, continuous and strictly increasing for t > 0, 

iic) it has an analytic one-to-one extension to an open set 11 C C with Qg <Z Q for some 
nonzero 6, and there is an a > such that lim(i^oo^ ^g^) t~q{t) = k, k > 0, 

iiic) Jq I In q(u)\du < 00 for all finite t > 0. 
With the above conditions we find. 

Lemma 3.7. Suppose ic)-iuc) hold then q^ e H{Qs) for e > and qo{t) = t«. Also 
for every sq > 0, q^ E C°°{^5 x (0,eo]), ^qe e C{Qs x [0,eo]) for i = 0,... and for 
every e G [0, eo] qe{t) is one-to-one on fls- qe{t) > 0, for (t, e) G (0, 00) x [0, eo] and 
I In q^(u)\du < 00 for all finite t > and e G [0, eo]. 
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Proof. Since q(t) e -ff (fij) we see that q(^) e C°° ((0, eo]) and from iic) that hnig^o Qit, e) — 
= uniformly on compact subsets of fij. This gives the continuity and analyticity 

properties of q{t,e) and shows that ^o(^) = . It follows from iic) that qe{t) > 0, for 
{t, e) G (0, cxo) X [0, eo]. That qe{t) is conformal for each e G [0, eo] and increasing for t > 
can be deduced from ic), iic), and the definition of qe{t)- To see the integrability condition 
write 



In q^{u)\du 



In 



e^/"q{u/e) 



eig(l/e) 

K ft/e 



du 

(3.26) 

/ , , / 1 . y s s , , / , , y 1 . y s ^ , , y , , 

< e 



/ \\n{s^^'^q{s))\ds + e \\n{s^ ' q{s))\ds + Ct{\\nt\ + 1) . 
Jo Jk 

Property iiic) implies the convergence of the first integral while iic) implies the integrability 

of the second term. □ 

From the definition of qe{t) we find that 

9r'(y) = er^(gQ + 0?/)-|- (3.27) 

Set = qf.{Qs)- The fact that q is conformal on Qs and maps (0,oo) into (0, oo) 
implies that = qeiins)+) C C+ and S'l = q,{{ns)-) C C_. 

Lemma 3.8. For each &xed e,q~^ G H{S'^). If y > r > then there is an eo such that 
q~^ > (j/, e) G [r, oo) x [0,eo]. If K is a compact subset of then there exists an e\ 
such that qj^ G H{K) for e G [0,ei], qj^ G C°^{K x (0,ei]) n C{K x [0,ei]). Finally 
Qe (^)x[ge(o),i](") _^ in L([0, 1]), where X[qe{o),i] tie characteristic function of the 

set [ge(0),l]. 

Proof. That q~^ G H{S^) is a conquence of Lemma 3.7. Since qQ^{r) = r'* > the 
continuity of q~^ in e implies that there is an eo so that q^^{r) > for all e G [0, eo]. The 
monotonicity of q~^ as a function of y shows that q~^{y) > for all {y, e) G [r, oo) x [0, eo]. 
Given K C Sq a, compact set let C O5 be such that K = qo{W). The continuity of go 
implies that W is compact. The uniform continuity of q^ (t) on compact subsets of x [0, eo] 
implies that there is an ei > and a compact set Wi C ^Is such that K C qe{Wi) for all 
e G [0, ei]. Thus q~^ is well defined on K for e G [0, ei]. The conformality of q^, Lemma 3.7, 
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and the inverse function theorem show that 57 e H{K)nC{Kx[0,ei])nC°°{Kx{0,ei]). 
To see the convergence we will show that the above family is uniformly integrable and 
then use Vitali's convergence theorem. Because of the continuity properties of this family 
only the uniform integrability near zero need be shown. Note that 

ln{q,{u))du^ lii{w)q;\wydw^tlii{q^{t))- ^ ' dw. 

Jo Jq,{Q) he{Q) 

Utilizing the fact that ^^(t) < 1 for < t < 1 we find, 

r-'''m<oM£H^rf„</Vaw)M«. 

Jo w Jq 

Lemma 3.7 shows that t^l°^ uniformly on compact subsets of R+ so that the 

uniform integrability follows from (3.26). □ 
A limit we will have use of is 

hm |y-«57i(y)|=6^fi + i) . (3.28) 

With the above assumptions we are able to give some representation formulas for Q. 
Suppose y > then for e sufficiently small q~^{y) > 0. Thus 



= 1, pL_ _ ^ + .(y _ AV^ 1 Lu 

^ ^ Jo \2aq^iu) 2a \l \2aq^{u) 2a J I 

= In f ^ - A + Jfj^-lY-l] q-\n:ydw. 

7g,(0) \2aw 2a \l \2aw 2a J j y J 

Integration by parts yields, 

'b+2a q-'^(w) dw 



(3.29) 



<5(2/,e) = y 



q~^(yu) du 



(3.30) 



'ifM u ^(1 - 6w)2 - (2aw)2 
Likewise for case 1 and y < we find 



Q{y,e) 



J a. 



'>-'^°- q^ ^{w) dw 
Qe(o) ^ Viy - - (2a«^)^ 
'^o-b q~^(—yu) du 



ImM. u ^(1 + 6w)2 - (2aw)2' 
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(3.31) 



The definition for V^{y, e) is 



' tJo \ 2aq,{u) 2a V \2aqe{u) 2a J ' 



(3.32) 



1 

H — / liiaqe{u)du, 
t Jo 

and the density for the equihbrium measure also simphfies to 

dp 



dy TTt J(y+)-i 



(7e+)-(.) ^i^+ip)-y){y-^r{p)) 
"^""^^ q~^{y/wydw 



(3.33) 



tn Jy(2a+b) u}^/4a'^ - (w - b)^ 



for cases 1-3 and max{(6 + 20)^^(0), 7^ (t)} < y < ^^{t). For case 1 and 7^ (t) < y < 
(6-2a)g,(0) 



dvf 1 /" dp 



dy Ji-rD-Hv) ^{^+{p)-y){y-^r{p)) 

1 q-H-y/wYdw 



(3.34) 



tTT J-y(b-2a) yj^J^a^ - (6 + i(;)2 

-It (*) 

while for case 3 and (6 + 2a)ge(0) < y < %{t) 

dul 1 f(^e)-\y) dp 



dy t^Ji.tr^iv) ^{^t{p)-y){y-i7{p)) 
1 q-^{y/w)'dw 



(3.35) 



A-2a t(;A/4a2 - (6 - tu)^ 



The above representations are obtained from Lemma 3.3 by setting u — q^ip), u — yv or 
—yv and then w — 1/v. 

Theorem 3.9. Suppose ic)-uic) hold and let r > 0. Then there exists an eo such that 

for all e < eo, Q{y,e) exists for y > r. Furthermore Q{y,e) has an analytic extension to 
S"^. If K <Z is a compact set, then there exists an ei such that Q{y,e) e H{K) for 
e e [0,ei], and Q e C°°{K x (0,ei]) n C{K x [0,ei]). For case 1 and y < -r there exists 
an eo such that for all e < eo, Q{y, e) exists for y < —r. Q{—y, e) has an analytic extension 
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to 5'^ and if K G is a compact set, then there exists an ei such that Q{—y, e) e H{K) 
for all e e [0, ei]. Also Q{-y, e) e C°°{K x (0, ei]) n C{K x [0, ei]). For each t e (0, tfi] 
and e > 0, V'{y,e) G if(C \ (-oo, 7+(t/„ e))), V^y^e) G C°°(C± x (0,t/,] x (0,eo]) 
for eo > 0. For each compact set K E C \ (— oo, 7*^(t/i, 0)) there is an ei such that 
tG eC°°{K X {0,tfi] X (0,ei]) and ^tV^ e C{K x (0,t/i] x [0,ei]), i > 0. Finaiiy tK[ 
and tVl have extensions to C+ and C_ respectiveJj tiat are in C{C± x (0, tfi] x [0, eo]) 
for eo > 0. 

Proof. From Lemma 3.8 we see that for y = r > there is an eo such that q^^{y) > 
for all e e [0,eo]. Since for fixed e, q^^{y) is an increasing function of y, q^^{y) > 
for all y > r. Thus Theorem 3.1 and Lemma 3.7 show that Q{y,e) is well defined for 
each e e [0, eo] and y > r and can be represented by (3.30) . Lemma 3.8 and Morera's 
Theorem, applied to the second integral in equation(3.30) gives an extension of Q{y, e) 
that is analytic on S"^ \ [0, q^{0){b + 2a)] for e e [0, eo]. Let K G he a, compact set. Then 
the above argument and Lemma 3.8 show that there is an ei such that Q{-, e) e H{K) for 
all e e [0, ei] and Q e C°°{K x (0, ei]). Using the convexity of S'^ chose ^e(O) < t so that 
the branch of the square root is respected and the line segment [t, ^q^] G 5^, e e [0, ei] so 
the second integral in (3.30) may be written as 

r's+2E q~^(yu) (ill fy rb+2E q~^[yu) du 

JniM u ^(l - buy - {2auY ~ JsiM J± u ^(1-6^)2 - (2aw)2 ' 

Lemma 3.8 now shows that Q G C(i^ x [0, ei]). An analogous argument applies for case 1 
and y < using (3.31). Since 7^(t) is an increasing function of t, Lemmas 3.2 and 3.7 
imply that V^{y,e) G H{C \ (— oo, 7"'"(t/i, e))) for e > and tV^ have extensions to C± 
that are in C{C± x (0,t/i] x [0, cq]) for eo > 0. Prom the definitions of ge(t), and and 
Lemma 3.7 we find that for each compact set K E C \ (— oo, 'y~^{tfi, 0)) there is an ei such 
that tV^ eC°°{K X {0,tfi\ X (0,ei]) and §^tV^ G C{K x {0,tfi] x [0,ei]) for i > 0. □ 
With the above results we are able consider extensions of the Langer transformation 
into the complex plane. We use the principal branch of the log. 
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Lemma 3.10. Suppose ic)-iiic) hold and b < 2a. For y > r > and t e (0,t/j], tfi < oo 
set 

Ci(t,y,e) = ^^^^^ g) _t)3/2 In (^2;(w,y,e) + V^(w,y,e)2 - Ij (3.36) 

witi t, e) = (t) ~ Tien there exists an sq such that Ci £ C((0, t/i] x [r, oo) x 
[0, eo]), is positive for (t, y, e) G (0, t/i] x [r, oo) x [0, eo] and ^Ci e C((0, t/i] x [r, oo) x [0, eo]) 
for aii i > 0. For £xed {t, e) e (0, t/j] x [0, eo], Ci(^) Vj ^) can be extended so that Ci{t, •, e) e 
H(S') with Ci e C((0,t/i] X S'). For fixed y e S% Re((t+(y,e) - t)3/2Ci(t, y, e)) is a 
decreasing function of t E (0,t/j]. For < e < eo there exists a 5i(e) > such that Ci 
is nonvanisiing for {t,y) E [tin, tfi] x Li{e) where Li{t) = {y : y E S'^ , |Imy| < If 
K (Z is compact there exists an such that ^ >^ K x [0,ex]), i > 0. 

Ifb<2a and —y<r>0 set, 

1 /-i^Cj/.e) 

C2{t,y,e) = —— ^ — / \n{z{u,y,e) + ^/z{u,y,e)^ -l)du. (3.37) 

(tp (i/,e) -ty/^ Jt 

with z(t, J/) e) = ^ — 2aq^{t) ' ^^^^ tiere exists an eo suci tiat ^2 £ ^"((0, t/i] x (— oo, — r] x 
[0,eo]), is positive for {t,y,e) E (0,t/i] x (-oo, -r] x [0, eo] and e C((0,t/i] x 

(— 00, — r] X [0,eo]) for aii z > 0. For fixed (t, e) e (0,t/i] x [0, eo], C,2{t,y,e) can he 
extended so that C2(t, ■,£) e H{-S') with C2 G C((0,t/i] x -S"). For fixed -y E S", 
Re(t~(y, e) — t)'^/^(^2(^) J/j e) is a decreasing function of t E (0,t/i]. For < e < eo 
tiere exists a 5i{e) > suci tiat Ci ^'s nonvanishing for {t,y) E [tin, tfi] x -£'2(e) wiere 
-£'2(e) = {y : y e — -S^, | Im j/| < 5i}. If K C —S^ is compact there exists an ek such that 
§,C2eC{{0,tfi]xKx[0,eK]),i>0. 

Proof. Set = In^ (^z{t,y,e) + ^ z{t,y, e)^ - 1^ with z{y,t,e) = - ^. For 

y > r > we see from Lemma 3.8 that there is an eo so that t^{y, e) > for aU (y, e) E 
[r,oo) X [0,eo]. Lemmas 2.1, 3.7, and 3.8 imply that p(t, y, e) = (j^t'^t^)^ e C'((0>*/i] x 
[r, 00) X [0, eo]) is positive in this region. From Lemmas 3.7 and 3.8 we see that 

00 

z{t,y,e) = l + ^Cj{y,e){t+{y,e) - ty , 
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where the sum is uniformly convergent for t in an interval about tp{y,e) with Cj{y,e) e 
C([r, cx)) X [0, eo]) nC°°([r, oo) x (0, eo]) for eo sufficiently small. Thus 2/, e) G 

C{{0,tfi] X [r, oo) X [0, eo]). For {t,y,e) G {0,tfi] x [r, oo) x [0, eo] set 

[ (t+(,,e)-t)3/^ j;i^'''\t+iy,e) -uy/Mu,y,e)du for t < t+(y,e) 
f{t,y,^)=\ ^p{t+{y,e),y,e) for t = t+{y,e) , 

[ ^t-t+{y,e))3/2 ^t+{y,e)(^ ~ ^P^y^ ^)y^Mu, y, e)du for t>t+{y,e) 

and for {t,y,e) e (0,t/j] x [r, oo) x [0, eo] set 

[ (4(,,,/_,)3/.+. j;^^'''\t+{y,e) -uy/^+%p{u,y,e)du fovt<t+{y,e) 
^/(^,2/,e) = < ^^§ip{t+{y,e),y,e) (or t = t+{y,e) ■ 

[ (,_4fa'e))3/2+i //+(y,6)(^ - e)y/^+%p{u,y,e)du for t > t+(y,e) 

Lemma 2.1, integration by parts and the mean value theorem for integrals shows that the 
definition of / is self consistent, / G C((0,t/i] x [r, oo) x [0, eo]), is positive in that region, 
and ^/ e C((0, t/i] x [r, oo) x [0, eo]). Now fix (t, e) such that 7+(t) > r and let y > 7+(t). 
We find from Theorem 3.1 and equation (3.32) that 

{t+{y, e) - tf/\i{t, y, e) = Q{y, e) - t{V\y, e) - k). (3.38) 

Since {b + 2a)q^{0) < (6 + 2a)q^{t) — 7"'~(t, e) we find from Theorem 3.9 that {tp{y, e) — 
t)^/^Ci(^) y, e) has an extension that is analytic in \ [0, 7"'"(t, e)] for each (t, e) G (0, t/^] x 
[0,eo]. Lemma 3.8 and the Schwarz reflection principal show that Ci{t,-,e) G H{S^) for 
{t, e) G (0, tfi) x [0, eo] and this combined with the above argument and Theorem 3.9 give the 
remaining continuity and differentiability properties of Ci . Taking the real part of the above 
equation then differentiating with respect to t yields — In \z{t, y, e) + \/ z(t, y, e)^ — 1| < 
which shows that for fixed y G 5''^, Re((t+(y, e) — t)^/^Ci) is a decreasing function of 
t. The dominated convergence theorem applied to the second integral in (3.30) shows 
that limj/^ooyes^ \y~'^Qiyi ^)\ > for < e < eo where equation (3.28) has been used. 
This coupled with (3.38) shows that there is a Si{e) such that Ci is nonvanishing on 
[tin.tfi] X Li(e). 

For b < 2a and y < if we replace above by = In^ (^z{t, y, e) + -^^(tTyT^P""^^ 
with z{t, J/) e) = ^ — 2aqe(t) arguments analogous to those above can be used to show that 
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C2 e C{{0,tfi] X (-00, -r] X [0,eo]), is positive for {t,y,e) e (0,t/j] x (-00, -r)] x [0, eo] 
and -§^(2 e C{{0,tfi] x (— cxd, — r] x [0, eo]). We now extend C2 to complex y. Fix {t, e) such 
that 77 (i) < —r and let y < j^it). Set G±{y, e) = G^^{y, e) — ItT ^^r. Then the extension 
of Lemma 3.2 to the case when e > (or by direct computation) we find that e) = 

G-{y, e) for y < 7"(t, e). Consequently there is a function e H{C \ [7~(t, e), 00)) such 
that = G± for y e C± and 

{t; {y, e) - tf/%{t, y, e) = Q{y, e) - tV*{y, e). (3.39) 

Theorem 3.9 now shows that {tp{y, e) — t)^/'^(2{t, y, e) has an extension that is analytic in 
—S^ for each (t, e) e (0, tfi] x [0, eo]. The rest of the proof follows as above. □ 
We now consider case 3. Again 7" becomes an obstacle to the smoothness of k^. 

Lemma 3.11. Suppose ic)-uic) hold and b > 2a. With y >r > and t e (0, tfi], let Ci be 
as in equation (3.36). For every {tfi, r), < tfi < (70")""*^ (''))) < r < 00, there exists an 
eo such that Ci £ C{{0, tfi] x [r, 00) x [0, cq]) is positive for (/;, e) G (0, t/^] x [r, cxo) x [0, cq] 
and ^Ci e C!{{0,tfi] x [r, 00) x [0, eo]) for all i > 0. For fixed (t, e) G (0,t/i] x [0,eo], 
Ci(t,y, e) can be extended so that Ci(t, •, e) G H{S^\[0,-f-{tfi)]) with Ci G G{{0,tfi] x {S^\ 
[0,7r(*/i)]))- -P^or iixed y e \ [0,%itfi)], Re{{t+{y,e) - tf/\i{t,y, e)) is a decreasing 
function oftE (0, t/^]. For < e < eo there exists a 5i{e) > such that Ci is nonvanishing 
for {t,y) G [t^n,tfi\ X Li(e) where Li(e) = {y : y g 5" \ [0, 77(t/,)], |Imy| < 5i}. If 
K G \ [0,7(7(t/i)] is compact there exists an ex such that (i G C((0,t/i] x K x [0, ex]) 
and £ C'llOj^/i] x K x [0, e^]), i > 0. Consider 

I rt 

C2{t, y, e) = - — — — — / \n{z{u, y, e) + y/z{u,y,ey - l)du (3.40) 

[t-tp {y,e)y/^ Jt-{y,e) 

with z{t,y,e) = ^ - 2aq,iu) ■ Let r > then for every [tin, tfi] x [^1,2/2] Vi < 2/2 with 
[z/i)2/2] G ['^j 00) ^J^d [Umtfi] £ ((7o')~"'^(z/2)) 00) tiere exists an eo suci tiat (2 is positive 

for {t,y,e) G [tin,tfi] x [2/1,2/2] x [0,eo], and ^(2 e C'([t,„,t/,] x [2/1,2/2] x [0,eo]) for 
aiJ i > 0. For fixed (t, e), C2 can be extended so that C2{trj^) €E H{{S'^ \ [7+(t),oo))) 
witi C2 e (^([tin,^/*] X {S^ \ [7+ (tin), 00))). -For eacL fixed y e \ [j+ (tin) , 00) , Re{t - 
tp (y, e))^/^C2(^5 2/) ^) an increasing function oftE [tin, tfi] . For < e < eo there exists a 
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S2{e) > such that (2 is nonvanishing on [tin,tfi] x I/2(e) where L2{e) = {y : y E , < 
Re(y) < 7^ (tin), I < S2}. If K C \ [joitin), 00) is compact there exists an ek 

such that -§^(2 e C{[tin,tfi] x K x [0, ex]), i > 0. 

Proof. The first part of the above Lemma follows as the first part of Lemma 3.10 noting 
that the restriction on tfi and \ [0, 77(t/i)] is necessary to avoid the singularity in 
due to 77. To show the second part of the Lemma fix r > then from Lemma 3.8 there 

exists an eo such that t~{y, e) exists and is greater than zero for all (y, e) G [r, 00) x [0, eo]- 
Set /c^ = ln'^{z{t,y,e) + \/ z{t, y, e)^ — 1) with z{t,y,e) = ^ ~ 2aq~{€)' th.3X 
t > t~{y, e) for all (t, y) e [tin, tfi] x [^1,2/2] in this case we increase j/2 as large as possible 
consistent with the restriction on tin- Lemmas 2.1, 3.7, and 3.8 imply that p{t,y,e) = 
(^JlShMilL^^ e C{[tin,tfi] X [2/1,2/2] X [0,eo]) is positive in this region. From Lemmas 3.7 
and 3.8, 2; has the expansion 

00 

z{t,y, e) = l + ^Cj{y,e){t-t-{y,e)y, 

where the sum is uniformly convergent for t in an interval about t~ {y, e) with Cj {y, e) e 
C([r, 00) X [0,eo]) n C°°([r, 00) x (0,eo]). This implies that §,p{t,y,e) e C{[tin,tfi] x 
[2/1,2/2] X [0, eo]) for aU i > 0. For (t, y, e) e [tin, tfi] x [2/1, 2/2] x [0, eo] and aU z > set 

[ (t-t-(y!e))3/2+i It-{y,e)('' " t' (y , e))^/^+' £,piu, y, e)du foTt>t-{y,e) 
^/(^,2/,e) = < 3^i§iP{tp{y,e),y,e) for t = t-{y,e) . 

[ (,-(,,e)-t)3m» J,'^^'^'\t;{y,e)-uy/'+^£,piu,y,e)du for t < t-{y,e) 

Lemma 2.1, integration by parts and the mean value theorem for integrals show that the 
definition of / is self consistent, / is positive for {t,y,e) e [tm,^/i] x [2/1,2/2] x [0, eo] and 
^/ G C{[tin,tfi] X [2/1,2/2] X [0, eo]) for all z > 0. To extend C2 to the complex plane fix 
{t, e). For r < y < 77 (t) set 

G±{y,e)^VUy,e)-ltTi7T. (3.41) 
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Then from Lemma 3.2 we find that G+{y., e) = G_ {y, e) for j/ < 7g (t). Consequently there 
is a function e H{C \ [7~(t, e), oo)) such that V* = G± for y e C± and 

{t - t-{y, e)f'H2{t, y, e) = tV\y, e) - Q{y, e). (3.42) 

Thus the extension of Lemma 3.2 to e > or direct computation shows that {t~{y, e) — 
t)'^^^C2(^i y, e) has an extension to 5''^ \ [77(^)5 00). Lemma 3.8 and the Schwarz reflection 
principle shows that C2 ^ H{S'^ \ [7+(t),oo)) for each (t, e) G [tin^^/i] x [0, eo]. The 
arguments above and Theorem 3.9 show that G C{[tin,tfi\ x (S''^ \ [7+(tin)5 co)) x 

[0, eo]) for all i > 0. Taking the real part of equation (3.42) then differentiating with 



respect to t yields ln\z{t,y,€) + \/ z{t^ y, e)^ — 1| > which shows that for fixed {y,e) e 
S"^ \ [7+(tiTj,, 00), Re((t~ — t)^/^C2) is an increasing function of t. The remaining part of 
the Lemma follows from the continuity and positivity of C2. □ 

Remark. We note that away from the turning point the above C, function are in fact 
C°°([tm,t/i] xKx (0,ei]) 

A consequence of the above lemmas is, 
Theorem 3.12. Suppose ii)-iiic) hold and b — 2a <0. For < Un < t < tfi set 

Pi{t,y,e) _ ,3^1,^ ^ 2/3 



, , -(^CHt,y,6))^/^ (3.43) 

Let L^{e) be as in the Lemmas 3.10 or 3.11. Then pi G H{Li{e)) for fixed {t, e) G [tin,tfi]x 
[0, ei] with pi G C((0, t/i] x Li{e)) . If K is a compact set in Li(0) then there is an ei such 
that G C([tin,i/i]xKx[0,ei]) for alii > 0. Also for {t,y, e) G [tin,tfi]x K x[0,ei], 

bp Z 

^ has an analytic extension so that, 

3 

Pi{t,y,e) = i-7vt di)t)i. (3.44) 
Likewise for b — 2a <0 and L'^{e) as in the Lemma 3.10 set 

f-f'':'\ =d<'{t,y,e)?'''. (3.45) 
tp [y, e)-t 2 
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Then p2 e H{L2{e)) for fixed {t,e) e [Un,tfi\ x [0, ei] with G C((0,t/i] x L^ie)) . If K 
is a compact set in L2(0) then there is an ei such that ^ttzt ^ C{[Un,tfi] x K x [0, ei]) 
for all i > 0. For {t, y, e) G [tin, tfi] x K x [0, ei], iias an analytic extension so that, 

3 f^^" "^^^ 2 

P2{t,y,e)^{-ntJ di)t)-^. (3.46) 

If 6 — 2a > tiieii for < tin < t < tfi -'et -^^^(e) be as in Lemma 3.11 and set 

P2{t,y,e) _ ^^^2u .^^2/3 



(^Cnt,2/,e)r • (3.47) 



t - tp {y, e) 2 

Tien p2 G H{L2{e)) for Bxed (t, e) G [tin,t/d x [0' ^i] with p2 G C([£in,t/i] x ^2(e). If 
is a compact set in L2(0) tiien there is an ei sucii tiiat G C([tiTi5 ^/i] x K x [0, ei]) 

for aii i > 0. For (t, y, e) G [tin; ^/i] x K x [0, ei], ^ iias an analytic extension so that, 

3 2 

P2{t,y,e) = {-i7r j dVl)-^. (3.48) 
In tiie above equations di>l and di>^ are given by Lemma 3.4. 



Proof. The continuity analyticity, and smoothness properties of pj, i = 1,2 foUow from 
the continuity, analyticity, smoothness, and nonvanishing properties of Ci, * = 1, 2 discussed 
in Lemmas 3.10 and 3.11 and the fact that Lj(0) i = 1,2 are simply connected regions in 
the plane. In order to show equation (3.44) note that for fixed t, and x > max{0,77(t)} 
equation (3.33) shows that ^ has an analytic extension to S"^ \ [0, 7^(t)]. Equation (3.44) 
now follows from (3.38) and the extension of Lemma 3.4 to e > . A similar argument 
for a; < follows using equation (3.34). Case 3 follows analogously from (3.42) and the 
extension of Lemma 3.4 to e > . □ 
For the extension of Theorem 2.6 to the complex plane we will use the solutions to 
the Airy difi'erential equation given by, Aio = Ai, and 

Ai±i(z) = Ai(^eT2^^/"^), (3.49) 

and the regions Sq = {z : \ argz] < f } and, 
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Since the above functions and Bi all satisfy the same differential equation we find, [35, 
p. 414], 

Ai(^e±2-/3) ^ le±-/3(Ai(^) ± iBi{z)). (3.50) 

Then it follows from the asymptotic expansions of Ai [35, p. 413] that Aij is recessive in 
Sj and dominant in S'j+i, and Sj-\ where the suflBx j is enumerated mod 3. Furthermore 
since the zeros of Ai are all real and negative [35, p. 418], Aii is non zero in U ^i. Two 
other solutions of the Airy equation that are non zero in the region of interest are tt;'^*^ 
which for complex values of z are defined by 

-^^Hy) = ^ (I) ' e^-'^'^^'-a^ (^y^e-) , i = 1, 2 (3.51) 

where we take the branch of the square root so that Re(2;^/^) > for 2; e and Re(2;^/^) < 
for 2; e which is the principal branch of 2;^/^. It follows from the properties of Hankel 
functions ([35, p. 238], [14]) that w^^^ is recessive in and dominant in S\ while tu^^^ is 
dominant in and furthermore neither vanish in S'o U S\ ([35], [14]) 
We rewrite g as 

p{i,y, e) 

a2 (t + |, e) sinh^(p5p')(t, e)/ 
and set, 

V'i(t, y, e) = ^Aio(e-2/3p(t, y, e)), (3.52) 

V'2(t, y, e) = ^Aii (e-2/3p(t, y, e)), (3.53) 

(t, y, e) = gw^^^ (e"'/ V(t, 2/, e)), w^^) ^ ^^ii (e-2/3p(t, y, e)), (3.54) 
^i(i)(t,y,6) = ei^(*'^'^)^«W(t,y,6), 

and 

«(2)(t, e) = e-i^(*'^'^)^ 

Set S' = S'o U 5*1 , and let O a region in the complex y plane. It follows from the asymptotic 
expansions of it* i = 1,2 [35] that it* are bounded functions. Furthermore both are non 
zero in U [14, 35, p. 254]. With the above notation we now have, 
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Theorem 3.13. Suppose (3.24) and ic)-iiic) hold. Let ai{t, e), bi{t, e) e C([0, oo) x [0, eo]) 
satisfy (2.24) with 6i(t, e) real and ai(t, e) strictly positive on every compact subset of 
(0, oo] X [0, eo]. In equations (3.52)-(3.54) let p = pi = {t+ - t)(Ci)^/^ and let Li(e) be 
given as in Lemmas 3.10 or 3.11 and L^{e) = Li(e)nC'+. Suppose K C L^(0), K compact, 
then there is an ex such that for each {y, e) e K x (0, ek] and all n : ne E [tin, tfi] there 
exists solutions fi i = 1,2 of equation (2.23) such that 

fi{n) = ipiin) + ri{n), (3.55) 

where 



riin) 




(n) 





fi{n) - i/jiin) 



< d{K)€, 1 = 1,2. (3.56) 



u(*)(n) 

Furthermore for fixed e, e H(Kr]C+) and E C{K x [0, e^]). Ifh-2a < let 

p = P2 = {t--t){C2ff^ equations (3.52)-(3.54). Ifb-2a>0 1etp = p2 = (t - t-)(C2)^/^- 
For both cases let L2{e) be given as in Lemmas 3.10 or 3.11 and i^J(e) = L2{e) fl . 
Suppose K C LJ(0), K compact, then there is an ex such that for each (y, e) E K x (0, ex] 
and all n : ne E [tin, tfi] there exists solutions fi i = 1,2 of equation (2.23) such that 

{-rfi{n) = Mn) + ri{n), (3.57) 

where 

, nin) ,_,(-im-)-^lM|<rf(K),, z = l,2. (3.58) 



•uW(n) ' wW(n) 
Furthermore for fixed e, E H{K n C+) and e C(K x [0, e^]). 

The proofs of the above Theorem closely follow the proof of Theorem 4.4 in [14]. 
However since the hypotheses are different a proof is sketched in Appendix A. 
We now have the important 

Lemma 3.14. With fi and /2 above, 

zes e 3 

ai((n + l)e, e) [/i(Wi)/2(tn) - /i(tn)/2(tn+i)] = 2^(1 + 0(e)). (3.59) 
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Proof. The Wronskian of Ai and Bi [35] is 

W^[Ai,Bi] = -. 

TT 

This coupled with the definition of Theorem 3.13, equation (3.50), and Lemma 3.5 in 
[14] (with u^^^ and tt^^^ replaced by /i and /2 respectively) give the result. □ 

IV. Singular Initial Value problem 

The previous results allow us to find solutions to the difference equation of a prescribed 
form that are valid uniformly in the neighborhood of a turning point. In order to obtain 
asymptotics for special functions we will need to find approximate solutions that satisfy 
the initial value problem. Unfortunately in the case when a{n) and \h{n)\ tend to infinity, 
in general the partial derivatives of a(t, e) with respect to t will not be bounded in a 
neighborhood of t = 0. This is because if a(t, e) and 6(t, e) are to be bounded then Ag 
must increase to infinity. (See section 6 for an example.) A similar problem arises in the 
case of varying recurrence coefficients [11], [25] which are connected to the continuum limit 
of the Toda lattice. While the main results. Theorem (4.4) below can be obtained from 
Theorem (5.5) of [14] the conditions on q allow simpler proofs which will be presented. We 
begin by switching to the equation satisfied by the polynomials pn = '2^Pn/kn, where k{n) 
the leading coefficient of Pn. Thus, 

Pn+l{x) + 2{bn - x)pn{x) + 4:alpn-l{x) = 0. (4.1) 

If the scaling indicated in the introduction is performed then writing p^ = Pn/ A" we arrive 
at the difference equation, 

Pn+i{y) + 2{b{ne, e) - y)pn{y) + ^o?{ne, e)pn-i{y) = 0. (4.2) 

The associated e-difference equation is, 

i;{t + e, y, e) + 2{b{t, e) - y)^{t, y, e) + 4a''{t, e)^{t - e, y, e) = 0. (4.3) 
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Following Deift and McLaughlin [11], Costin and Costin [7], and [14] we look for approxi- 
mate solutions of (1.2) of the form 



If we expand ilj{t±e) in powers of e then equate like powers the eikonal (e°) equation gives 
(see section 6), 

e-oW + 2{b{t, e)-y) + 4a{t, efe'^'o^^^ = 0, (4.4) 
while the geometrical optics equation (e^) gives 

Combining 

(e*°(*))' - 4:a{tfe-'°^'yy = {e'^^'^ + 4:a{tfe-'°^'y)s'^{t) - 8a{t)a'{t)e-'°^'\ 

with the derivative of (4.4) yields, 

(e-°(^--)'-4a(t,6)^e-o(^--)^)- b{t,ey so{t,e)" 

' ^ 2(e*o(*'^)' -a(t,e)2e-*o(*'^)') e^o(*'^)') -a(t,e)2e-^o(*'^)' 2 ' 

Solutions to the above equations are given by, 

4 {t, e) = y In - b{u, e) ± {y - b{u, e)f - ^f) du, (4.7) 

and 

4(t,e) = -^ln((y-6(t,e))2-4a(te)2) 



+ i ln(y - b{t, e) ± ,/{y - b{t,e)r - 4a{t,er) 
1 b'{u, e)du 



(4.8) 



^27 ^{y-b{u,t)Y-Aa{u,ty 
A formula which we will have use of later is, 



±.,y,_ h'{t) 2a\t)' 
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where 

hi{t, y, e) = v/(2/-K^,e))2-4a(t,e)2, (4.10) 

and 

hf{t, y,e) = y- b{t, e) ± V(y - K^, e))^ - 4a(t, e)^ (4.11) 
We will take as our approximate solutions 

, , , , h+(t,y,ey/^ ( b'(u,e)du \ 

hi[t,y,ey/^ \Jo 2hi{u,y,e)J ^^^^^ 

X exp ^1/e In /ij(u, J/, e)(iu^ , 

V'l y, e) = -TTT X exp - / — — 

X exp I 1/e / ln/i^(it, y,e)(iit I . 

\ -/tin / 



and, 



Let, 



(4.13) 



7±(t, e) = 6(t, e) ±2a(t, e), 

A{i, e) = inf 7~(t, e), B(t, e) = sup7''~(t, e), 
[o,t] [o,t] 

e) = (-00, B{i, e)], and /(f, e) = [^(t, e), B{i, e)]. 

We now examine the analytic properties of the above functions. 

Lemma 4.1. Suppose a{t,e) and b{t,e) are given by equation (3.24) and conditions ic) 
and iic) hold. Then for each (t, e) e [0, tfi] x [0, eo], In/iJ is nonvanishing and analytic on 
C\ (— oo, ^~^{t, e)] while and hi are nonvanishing and analytic onC\ [7~(t, e), ^~^{t, e)]. 
For each (t, e) e [t^n, ^/i] x [0, eo], Inh^ is analytic and nonvanishing onC \ (— oo, 7"''(t, e)] 
and h2 is analytic and nonvanishing on C\ [7" {t, e), ^'^{t, e)] . For any compact set K C C\ 
Isitfi^ 0) there is an ex such that In /ij G C([0, t/^] x i^T x [0, e^^]) and In G C([^m5 tfi] x 
iiT X [0,6^]). For any compact set K C C \ I{tfi^Q) there is an such that /ij, hi E 
C([0,t/i] X K X [0,e^]) and §,hi, ^h+, ^lnh+ e C{{0,tfi] xK x [0,e^]), i > 0. 
Likewise there is an such that ■§^h2 , and ^t+t In/i^ e (^([t^n, ^/i] x K x [0, e^]), i > 0. 
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Proof. The continuity, differentiability and analyticity properties of In h2 follow from 
equation (3.13), the fact that ln(^ — e [7~(t, e), 7"'"(t, e)] is analytic for 2; e C \ 

{—00, 7"'~(t, e)) and the differentiability part of Lemma 3.7. The continuity, differentiability, 
analyticity and nonvanishing properties of /i^ and hi follow from their definitions (see 
equations (2.8) and (2.9)). The fact that a{t,e) > for t > and the formula h2 = 
imply the results for this function and its log. □ 
In order to simplify the analysis below we make the assumptions that there exist 
constants Ci > and C2 > such that, 

ivc) \§jq{ty\ <Ci{t + C2)P^''^\ j = 0, . . .2, i = 1, 2, t > where 



a •' a •' ' 

ii±-j = -i 



Set 



-j^ /-tie /-tie 

J%+{t) = — s^{uy"{t±e- ufdu+ s^{u)"{t±€-u)du, 
«o 2e Jt Jt 



and 



-j^ /-tie /-tie 

+{t) = - / s^{u)"{t±e-u)du+ / s'liuydu. 

Lemma 4.2. [Lemma 5.2, 14 errata] Suppose a{t, e) and b{t, e) are given by equation (3.24), 
conditions ic), iic) and ivc) hold, and < eo <^ tfi. Then for < e < eo and (t, y) G 
[0, t/j] X C \ lB{tfi,e), is nonzero and for fixed {y, e) satisfies the difference equation, 

^+it + e) + 2{h{t,e)-y)^t{t) 

(4.14) 

+ Aa^{t, e)tl;t{t - e) = ?7+(t)V'^ (^)- e<t<tfi-e 

For eo sufRciently small, ip^ e C{[0,tfi] x (0,eo]) x H{C\ Isitfi, e)) and rj^ e C{W+) x 
H(C \ iBitfi, e)), W+ = {{t, e) : e e [0, eo], e <t < tfi - e}. For each compact set K C 
C\lB{tfi, 0) there exists an sk such that i;^ G C([0, tfi] x K x (0, ex]), Vi e C{W+^ x K) 
and 

J2 |?7+(ie,y,e)| < cmax(6e^/",e2/«,e) (4.15). 

i:iee[e,tfi] 

Here c depends on K and tfi. If t = tn = ne then the sets Isitfi, e) and Isitfi, 0) in the 
above statements maybe replaced by the sets I{tfi, e) and /(t/i, 0) respectively. Likewise 
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for each < e < eo and (t, y) e \tinitfi\ x C\ lB{tfi, e), '4>i is nonzero and for fixed (y, e) 
satisfies the above difference equation for tin + e < t < tfi — e with ry^ replaced by, r]^ . 

G C([t,„,t/,] X (0,eo]) X H{C\lB{tn,e)) and 77+ G C{W-) x if(C \ /^(t/,, e)), W~ = 
{{t, e) : e G [0, eo], tin + ^ t tfi — e}- For each compact set K C C \ Isitfi, 0) there 
exists an ck such that ip^ G C{[tin, tfi] x K x (0, ex]), rji ^ ^(^e^ ^ 

X] l^r(^e,y,e)| < ce, (4.16) 
witi c depending on K and the interval [tin, tfi]. 

Proof. The conditions on q imply sup^gjo.eo] /o^' < Thus Morera's theorem, 

equations (4.10), (4.11) and Lemma 4.1 show that i/jf exist and have the claimed continuity 
and analyticity properties. Using Taylor's theorem with remainder, equations (4.4) and 
(4.6) we find (suppressing all variables but t). 



r;+(t) = e^^W'(j:+,.W + i?:.,+ (^)) 



+ Aa{t, efe-^o W (J- (t) + R;^ (t)). 



(4.17) 



where 

i?± ^+(t)) = L± ^{tf [ exp(L± ^At)u){l - u)du. 



The properties of 77^*" except equation (4.15) follow from Lemmas 4.1 and condition ivc). 
If t = tn = ne then for y < A{tfi, e) we find 

- \nh^{t,y±,e)dt= - I \n\h'^{t,y,e)\dt±in-K. (4.18) 
^ Jq ^ Jo 

Thus 'ilji{tn, y+, e) = ipi{tn, y—, e) which along with Lemma 4.1 gives the stated properties 
with Isitfii e) and Isitfi, 0) above replaced by I{tfi, e) and /(t/i, 0) respectively. To obtain 
(4.15) note that 
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which when combined with the definition of J _,_ +, i? + + and Lemma 4.1 give (4.15). 
Since 



Vi 



it) = e^o-W'(jt ^-{t) + Rt- .-(t)) +4a(t,6)V^o-W'(j- (t) + R- (t)). (4.19) 



Sq ,s 



an argument similar to the one above gives the result for i/;^ and rj^ . Equation (4.16) 
follows since tin > 0. □ 
As we will use perturbation theory we suppose for each e e [0, eo], 

sup |6(t,e) - 6i(t,e)| = O(e^) = sup \a'^{t,e) - al{t,e)\. (4.20) 
te[o,tfi] t£[o,tfi] 

With the above we obtain solutions to the initial value problem. 

Theorem 4.3. [14, Theorem 5.5] Suppose a{t, e) and b{t, e) are given by equation (3.24), 
and ic), iic), and ivc) hold. Let ai{t,e), bi{t,e) G C([0, cxd) x [0, eo]) satisfy (4.20) with 
bi{t,e) real and ai{t,e) strictly positive on every compact subset of (0, oo) x [0, eo]. Let 
y e C\I{tfi, 0) and Pn{y, e) be a solution of 

Pn+i(2/, e) + 2(6i(ne,e) - y)pn{y, e) + 4ai{ne, efpn-i{y,e) = 0, 

(4.21) 

Po(y,e) = 1, Pi(y,e) = 2(y-6i(0,e)). 

Then there exists an ei such that j/ e C \ /(t/i, e), < e < ei. For each e, < e < ei and 
all n : ne & [0, tfi] 

Pn{y,e) ^i;+{n,y,e)/i;+{0,y,e){l + Mn,y,e)), (4.22) 

with tjj^ given by equation (4.12) and 

|0i(n,y,e)| < dmax(6e^/", ae^/", e). (4.23) 

If K is a compact set in C \ /(t/i, 0) and y & K then there is an €k such that (4.23) holds 
uniformly on K x [0, sk] and 0i(n) e H{K) n C{K x [0, sk])- 
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Proof. From condition ivc) and equation (4.20) we find that a(0,e)^ < ce«, \bi{ie,e) — 
b{ie,e)\ < c\e^ and |af(ze, e) — a^(ie, e)| < C2e^. This coupled with equations (4.15) and 
(4.16) and equation (5.28) of [14] give equation (4.23) (note that in equation (5.28) L(l, n — 
1) should be replaced by S{n — 1), see [15]). The smoothness properties for 0i follow from 
Lemma 4.2 and equation (5.30) in [14] with the correspondence /+ ipf. 

V. Matching 

The previous results can be combined to obtain a uniform asymptotics for solutions of 
the initial value problem. This is accomplished by matching solutions in various overlap- 
ping regions. Throughout this section we will suppose ai(t, e), 6i(t, e) G C([0, oo) x [0, eo]) 
with bi{t, e) real and ai{t, e) strictly positive on (0, T) x [0, eo] T > 1. 

Lemma 5.1. Suppose that conditions iic) and ivc) hold then 



where 



and 



|Cn|=0(-), 

n 



n = £^ln (1^) du - j^{B2 - B^iu - [u]))-^lnq{u)du (5.1) 

Here B2 and B2{x) are the second BernouUi number and polynomial respectively. Also if 
equation (2.24) holds then for n such that tin < tn < tfi 

«iH = n -^7^ = /^i(l + 0(e)), (5.2) 

where 

-1 = n 4^- 

f=i ai(ze,e) 

Proof. Conditions ic) and ivc) show that the second integral in k converges. The result 
follows by using Euler-Maclaurin formula [35, p. 285]. Prom (2.24) we see that for ti > 
tin, In -^^^ =0(€^) which implies (5.2). □ 

""1 i^J 

With the above we now prove, 
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Theorem 5.2. Suppose a{t, e) and b{t, e) are given by equation (3.24) and conditions ic)- 
ivc) hold. Let ai(t, e), e) satisfy equations (2.24) and (4.20). Suppose Pn{y , e) is such 
that pn satisfies the initial value problem (4.21). ffyE C+ \ {0} then there exists an ei 
and ats ^ (0) V«] ^^ch that for each < e < ei and all n : ne E [0, t^], 

Pn{y, e) = 

( (,-K^„..)/-4a^(.„ + j..) ) * - ■"<"<-'-^/^^^'- (5.4) 

X (1 + ^i(n,y,e)), 

wiere K(n) and Ki(n) are given in Lemma 5.1, t„ = ne, zi(u, j/, e) = '2a(u^&l^ ; 

12 

|0i(n, y, e)| < (imax(|6|e« , e« , e). (5-5) 

If K is a compact set of C+ \ {0} witii y E K then there exists an so tiiat tiie above 
error can be made uniform on K x [0, e^] for aii n : ne E [0,tB]- If V & ^i^(O) given in 
Lemmas 3.10 or 3.11 then there exists a tin < ts and an €2 < ei such that for < e < €2 

and all n : ne E [tin, ifi] 

^iQi^y^l = ««i2y7ry^i(ne,2/,e)(Ai(e"^pi(ne,2/,e)) + ri(n,y,e)), (5.6) 

where Q{y, e) is given in equation (3.30), k by (5.1), ki by (5.3), and pi by equation (3.43), 
and ^ 

- u \ ( e"^Pi(^n,y,e) V . 

'"''"'^''^=i((!/-6(t„..))^-4a(«„+|,.«j ■ («-^) 

The error term satisfies 

|ri(n,2/,e)e^^^(^?(*"'^'^»| < dmax(|6|ei , e« , e). (5.8) 



If K is a compact set in i^^(O) tiien there exists an €k such that the above bound is 
uniform on K x [0, €k] for all n : ne E [tin, tfi]- 

If b — 2a < and y E i^J(O) given by Lemma 3.10, then there exists a tin < ts 
an 62 < ei such that for < e < 62 and all n : ne E [tin, tfi] 



Pn{y,e) 
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= {-l)'^KKl2^/7^{^g2{ne,y,e) (^Ai(e 3p2(ne, y, e)) + r2(n, y, e)) , (5.9) 



where Q{y,e) is given in equation (3.31) and 



^^''"■''■''=^(fe-6(t„,e))^-4a((„+|,enl ' 

The error term satisfies 

|r2(n,2/,e)e^^*^(^2^(*")'^'^)| < dmax(|6|ei , , e). (5.11) 



If K is a compact set in i^J(O) then there exists an ex such that the above bound is 
uniform on K x [0, e^] for all n : ne E [tin, tfi]- 

Proof. See Figure 2 for a visualization of the above times. Fix y in C+ \ {0}. Since 
7o"(0) = we see from the definition of I{t, e) that there exists a is £ (0, tfi] such that for 
all (t, e) e [0, ts], y E C+ \ I{tB, 0). Thus Theorem 4.3 says there exits an ei so that for 
each e, < e < ei and all n : ne e [0, t^], 

Pn{y, e) = ^^tiu,v,e)'+sUu,y,ey)du^^ ^ ^^^^ ^^^^ ^5_^2) 

where tn = ne and equations (4.7)-(4.12) have been used. We find using Taylor's Theorem 
that, 



e / ~ ^(2/ - i>(w, e))2 -Aa[u + |, j 

1 ^, . ,^ 1 /•*" la(»,f)a'(»,r) , 1 , 

= -/ ^T^hl{u,y,e)du- - —r- — — -du+- ei{u,y,e)du, 

e Jo 2 Jq h^{u,y,e)hi{u,y,e) e Jq 



(5.13) 



and 



ln((y - h{tn, e)f - Aa'^{tn + |, e)) = ln/ii(t„, e) + e2(tn, y, e)- (5.14) 



Condition ivc) shows that, 

|ei| < c(y) max(6e"'^/'*, ae« , e), (5.15) 
and condition iic) shows that for strictly greater than zero. 



e2(tn,y,e)| < c{y)e. 
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(5.16) 



It also follows from conditions ivc) that 



/ii(0,y,e) = y(l + e3(0,y,e)), (5.17) 

with 

1 2 

|e3(0,|/, e)| < c(|/) max(6e« , ae« ). (5.18) 

Substituting (4.9) into equation (4.6) for sf' then using equations (5.13), (5.14), and (5.17) 
in (5.12) yields. 



{y-b{tn,e))^ -4:a{tn + ^,e)^J (5.19) 
X C HM^,y,e)+VMu,y,e)--i)du^^ ^ ^^^^^ ^^^^ 

where 0i has the same properties as and thus satisfies the bound (5.5) by Theorem 4.3. 
This yields equation (5.4). If y G -£'^"(0) then Theorem 3.12 shows that there exist a 
Un < ts and an €2 so that y e Lf^e) for all < e < ei and pi is well defined for 
(t, e) G [tin.tfi] X [0, 62]. For fixed e the function on the right hand side of equation (5.4) 
satisfies the same difi'erence equation (2.23) as Pn{y,^) so that for all n : ne E [tm,^/i] 
Theorem 3.13 implies that it can be written as ci/i + 02/2- In order to compute these 
coefficients restrict n : ne & [tin, ts] in which case equation (5.4) can be used for p^- For n 
in this region Lemma 5.1 shows that K,{n) = k(1 + 0(e)) , Ki(n) = ki(1 + 0(e)) and from 
equations (3.38) and (3.43) we find that 



ei<JM ' V(!/-i'(t«,e))'-Mtn + f,e)^ 

Since there is no turning point in [tin^^s] we may use the asymptotic expansion of the 
Airy functions. 



1 



M{e--^p{tn)) = ^ e-s-Mtr.)- (1 + 0(6)), 



and 

€6 



Bi(6-^p(y ) = , e^7P(^-)^ (1 + 0(6)), 
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in /i and /2. Using the Casorati determinant Ca, 



Ca[tti, U2]{n) = ai{ne, e){ui{n)u2{n — 1) — ui{n — l)u2{n)), (5.20) 

we find '^^^ca^/i'^alCnf and ^'^Sf/^tlKn)' '^"^ ' The above equations coupled with equa- 
tions (3.56) and (3.59), show that 

ci = y \ \^r\nm,y,e)\ (5.21) 

and, 

where Um is the smallest value of n such that ne e [tm^^B] and the fact that Ca is 
independent of n has been used to obtain the above equations. Theorems 4.3 and 3.13 
show that 

1 2 

|f*| < c(2/) max(|6|e« , e~, e), i = l,2 
Since the real part of pi(t)^/^ is a decreasing function of t we find, 

IC2/2I < c(y)e-i^«-«(^?''«"-+i)^'^'^))|f2(n^,y,6)|, (5.22) 

which is exponentially small. This gives equation (5.6). The uniformity of the error terms 
follows from Theorems 3.13 and 4.3. 

If 6 - 2a < and ?/ e L+(0) where L2(0) is given by Lemma 3.10, there exists an €2 
such that y G -^2^(e) for all e G [0, £2] so that p2 is well defined for (t, e) G [tin, ^/i] x [0, €2]. 
For t„ G [0,tB] we find from equation (4.18), 

Pn(y,e) = 

X (1 + (/)2(n,2/,e)), 

(5.23) 

Prom Lemma 5.1, equations (3.39) and (3.45) we find for fixed (y, e) and all n : ne G [tin, t/i] 



^ ^ ^ xe"3i^2(t„)2 ^;L + 02(n,u,e)), 

ehQ^y^-) ' V((y-K*n,6))'-4a(tn+f,6)2,e),' ^ 



1 



(-i)>n(y,6) / 
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where Q{y, e) is given by equation (3.31). In this region we also have from the second part 
on Theorem 3.13 that 

e^Qiv,^) - ^Mn) + C2f2{n) 
where = g2{tn){Ao{e-^/^P2{tn)) + r^t^)) and /2(t„) = g2{tn){Ai{e-^/''p2{tn)) + 

i^'^{tn)) with g2 described in equation (5.10). Since there is no turning point in [tin, ^a] the 
asymptotic expansion of the Airy functions can be used so that the coefficients ci and C2 
may now be computed using the Casorati determinant as above. The uniformity of the 
error terms follows from the second part of Theorem 3.13 and Theorem 4.3. □ 







y 














^^^^^ t;(y) tfi ' 



X y 

Figure 2: Various times used in Theorem 5.2 
A case where the above results can be applied is when 6 = and let T2i{x) be a 
monic even polynomial in x with real coefficients and let T2i{xo) = and T2i{x) > for 
all X > xq. Let q{t) = T2i^{t) be the inverse branch of T21 that is positive for all t large 
enough. From the assumptions on T21 we see that q{t) is monotonically increasing for i > 
and since x = 00 is a critical point of order 21 the open mapping theorem (see Rudin[38] 

p), shows that in a neighborhood of infinity q{t) has the piuseiix series representation 

00 

n=-l 
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Furthermore q{t) has an extension to a wedge {t : argt < 5} where 5 is the argument of 
the closest non zero critical value to the positive real axis. Because of the conditions on 
T21 and T21 there are no critical values on the positive real axis. If t = is a non analytic 
point of q{t) the the piuseux series about t = of form Yl'^=i Cnt^^^ with p > 21. From 
the above argument we see that. 

Lemma 5.3. Suppose q{t) is the inverse branch of an even monic polynomial with the 
properties described above. Then q{t) will satisfy ic-ivc) with a = 21. 

We now consider the case when b — 2a > and < y < 7^(t/i) then we will be in 
the two turning point case and besides ts above we will need some other times in order 
to separate these points (see Figure 3). We will restrict ourselves to the real line since at 
this time we are unable to control the error terms for the matching problem beyond the 
first turning point in the complex plane. 

Theorem 5.4. Suppose a{t, e) and b{t, e) are given by equation (3.24) with b — 2a > 
and conditions ic)-ivc) hold. Let ai(t, e), 6i(t, e) satisfy equations (2.24) and (4.20) and 
< y < ^o{tfi). Then there exist Ib-, tin, tin, tfi and an e\ satisfying the inequalities 
< tin < tB < tin < tfi < ^fi ^ucli that t+{y,e) ^ [0,tB], t+{y,e) e {tin,ifi), t-{y,e) ^ 
[tin,ifi], t~{y,e) e {tin, tfi) and t+{y,e) ^ [iin,tfi] for all < e < ei. Suppose Pn{y,e) is 
such that pn satisfies the initial value problem (4.21). Then for < e < ei, ei sufEciently 
small and all n : ne e [0, t^], Pn{y,^) satisfies 

Pniy,e) = 

^ ' ' \{y-b{tn,e)Y -Aa^{tn+'^,e)) 
X (1 + 0i(n,2/,e)), 

where K,(n) and Ki(n) are given in Lemma 5.1, tn = ne, zi(u, y, e) = 2a{u+^% ' 



12 

')i{n,y,e)\ < (imax(|6|e« , e« , e). (5.26) 



For all n : ne & [tin,tfi 



= i^iii2V^9i{ne,y,e){Ai{e-^ p^{n€,y,€)) + ri{n,y,€)), (5.27) 
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where Q{y,e) is given in equation (3.30), k by (5.1), pi by equation (3.43), gi by equa- 
tion (5.7), and the error term satisfies 



3 

(n,|/,e)e*^^(^?(*"'2''^))| < dmax(|6|e^ , e« , e). (5.28) 



ri 



For all n : ne & [tin,tfi] 
Pn{y,e) 

= (-l)'^2KKi0ry^2(ne,y, e) 



smQr,(y)^ +Ri{y,€)^ (Ai(e-2/3p2(t„,y,e)) 
+ fi(tn,2/,e)) +2 (^cos Qr,(2/)^ +i?2(2/,e)^ Bi(e-2/3p2(tn, 2/, e)) + r2(tn, 2/, e)) 



where 
Ifere 

|fi(ne,y,e)eA^^(^l(*"'^'^))|, |f2(ne, y, e)e-^ ^^(^^*-^'^))|, |i?i(|/, e)|, ||i?2(2/, e)| 

1 2 

< cimax(|6|e« , e« , e). 

The error terms in equation (5.28) and (5.30) are uniform on closed intervals of [0, oo). 

Proof. See Figure 3 for a visualization of above times. Equations (5.25) and (5.26) follow 
from equations (5.4) and (5.5) of Theorem 5.2. Likewise if in Lemma 3.11 we replace tfi in 
Li(e) with tfi, equations (5.27) and (5.28) also follow from Theorem 5.2. We now consider 
the interval [tin, ^/i]- In this interval we have by Theorems 2.6 and 3.13 that 

with /i = (72(tn)(Ai(e-2/3p2(tn)) + ri(tn)) and ^(tn) = (72(tn)(Bi(e-2/3p2(tn)) + r2(tn)). 

By selection the interval [tin, tfi] is void of turning points so for small enough ei the 
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asymptotics of the Airy functions on the negative real axis, may be used [35, p. 392]. Thus 
with n : ne e [tin,tfi] equation (5.27) can be recast as 

Pn{y,e) 



^Q{y,e) 



= 2««i ( o /, . e n V ) ( ( |-(-Pi(^e>y>e))' - 7 ) +riin,y,e)] . 

\2a{tn + ^,e)smzi{tn,y,e)J \ \3e 4 J ) 

(5.32) 

Equations (2.1) and (2.17) show that 

^(-Pi(tn))^ = cos-1 % d« + 7r(tn - V(|y)) + -(-p2(t,)) 2 . 

Using equations (3.24), (3.25) and relations between and , we find 



y 

1 1 y- K'^) \ J -1/ M -I ( y-bw 



cos ( ——, -T- ] du = (to)' cos ( — 1 dw 

4iy) V2a(«+i); ^ ^ V 2a«; ' 

b+2a 



du TT 



V4a2 - (w - 6)2 2' 

where integration by parts and (3.27) have been used to obtain the last equality in the 
above equation. Since t„ = ne we find 

cos (J- {-p,{ne,y,e))i - ^) = (-l)-sin (-p2(tn))^ + Jr,(y) - ^) . (5.33) 

where Fg is given by equation (5.29). Therefore 

= 2««i ( ^ ^ e ^- n \Y (sin(|-(-P2(^n))n-r,(y)-^) +fi(n,y,e)) 

V2a(tn + f,e)sin2;i(tn,y,e)/ V V^e e 4/ / 

(5.34) 

Likewise in this interval 

h{tn) = I ^ u • ' (, ^^m f|-(-P2(tn))'/' - 7) +ri(tn)) , 

(7r2a(tn+ f,e)smz2(tn,2/,e))^/^ V V3e 4/ / 



and 



C'/^* /. /2, , 3/2 TT 



f2{tn) = -^TTVTT^^^^ ^ TTTTi sin (-p2(t„))3/2 _ + f2(tn) • 

2(7r2a(t„+ |,e)sin2;2(tn,y,e))^/^ V V3e 4' ' 
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The error term fi(t„) satisfies equation (5.11) and fi(t^) and f2(tn) satisfy a similar bound 

3 3 

with Re(pf ) replaced by ^R-eCpl) respectively. 

Comparison of the coefficients in equations (5.31) and (5.34) (or using the Casorati 
determinant) for ci and C2 and utilizing the fact that sin^i = smz2 yields 



and 



C2 = -4kki^^^ I^cos Qr^d/)^ + i?2^ , 



where |-Ri| < (imax(6ea , €« , e), i = 1, 2 which follows since we are in the oscillatory region. 
The result now follows. □ 




Figure 3: Various times used in Theorem 5.4 
Of special interest is when e = 1/N and t„ = 1 which implies that n = N. From 
equation (3.25) we see that q'e(l) = 1 and 7^ = 6 ± 2a. The turning points are at values 
of y such that y — b + 2aory — b — 2a. In the results below we suppress the dependence 
upon t since t = 1. 

Theorem 5.5. Suppose that b — 2a < and t = 1. For y e C+ \ [6 — 2a, 6 + 2a] 

PN{y,l/N) = 

KKi 



(y - by — Aa"^ 
The bound 



Jo H^ii^,yA/N) + ^zi{u,y,l/N)^-l)du^-^ ^ ^^^^^ 1/N)). 



\My, ^/N)\ < max(VA^i/«, l/AT^/", 1/N), 



(5.36) 
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holds on compact subsets of C+ \[b — 2a, 6 + 2a]. For y e -£'i'(0) given by Lemma 3.10 and 
N sufRciently large, 

PN{yA/N) ^ 

f,NQiy,l/N) 

1/4 3T) 

where Q is given by equation (3.30), pi by (3.43), n by (5.1) and ki by (5.3). The error 
term ri satisGes the bound 



|ri(2/,l/iV)e3^^^(^i^^'^/^))| <dmax(|6|e-,e-,e), (5.38) 

uniformly on compact subsets of Lf{0) . Ifb — 2a<0 and y e I' J (0) given by Lemma 3.10 
then for N sufEciently large, 

(-l)>Ar(y,l/iV) _ 

giVQ(j/,l/iV) 

2/3 1/2 (5-39) 

2^^i^AK^ (Ai(iV2/3p2(2/, 1/iV)) + r,{y, 1/N)), 

where Q is given by equation (3.31) and p2 by (3.45). The error term r2 satisfies the bound 
|r2(2/,l/iV)e^^^^(^2 (^'^/^»| <dmax(|6|ei,e*,e), (5.40) 
uniformly on compact subsets of L^(0). 

Proof. For y e (C+ \ [6 — 2a, 6 + 2a] as t„ varies from to 1 no turning points are 
encountered. Thus the result follows from Lemma 5.1, equation (5.25) of Theorem 5.2 
with tn = 1 and the analytic properties of the solution. The uniformity of the error 
also follows from Theorem 5.2. The result for the remaining regions follows in a similar 
manner. □ 
Likewise for 6 — 2a > from Theorem 5.4 with = 1 we find. 
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Theorem 5.6. Suppose that 6 - 2a > and t = l. For y e C+\[0,b + 2a] 
PN{y,l/N) = 



( Y /^'ln(«i(u,y,l/iV)+V«lK2/,l/^)'-l)d"/;L + (/)i(|/,iV)), 

V(y-6)2-4a2y ^ 



where (j)i satisGes equation (5.36) uniformly on compact subsets of C+ \ [0, 6 + 2a]. For 
y e -^i^(O) given by Lemma 3.11 and N sufficiently large, 

PN{y, 1/N) ^ 

^NQ{y,l/N) ^ ^ 

'^""'^ ( ^-^^-iy ) (Ai(iV^/^pi(y, 1/N)) + T,{y, 1/iV)), 

where Q is given by equation (3.30), pi by (3.43), n by (5.1) and (5.3) by (5.3). The error 
term r\ satisfies (5.38) uniformly on compact subsets of Lf{0). For ye (0, 6 + 2a) 

{~l)^PN{y,l/N) 

^NQ{y,l/N) 

^ ^[ NV^p^{yMN) \''\^ 
= '""^^ Uy-6)^-4a^ ) ^ (5.43) 

"(sin(iVri/^(y)) + i?i(y, l/iV))(Ai(iV2/3p2(y, l/AT)) + fi(y, 1/iV)) 
+ 2(cos(Arri/^(y)) + i?2(y, l/iV))(Bi(Ar2/3p2(2/, at)) + f2(y, 1/N)) 

where uniformly on compact subintervals of (0, h + 2a) 

|fi(y,l/iV)ei^^^(''?(^'i/^))|, |f2(y,l/iV)e-i^^^(''2^(^'i/^))|, 
|i?i(y, l/iV)|, |i?2(y, l/iV)| < cZmax(|6|/iV^, 1/iV*, 1/iV). 



(5.44) 



Remark. In the language of Deift et. al [10] and Biak et. al [2] the regions {h + 2a, oo), 
(6 — 2a, h + 2a) and (0, h — 2a) are caUed bulk, band and saturated region respectively. 
The neighborhoods around b + 2a are called a band/ void edge while the neighborhood of 
6 — 2a is called a band/saturated edge. Due to the error of the matching problem in the 
above technique the exponentially small error obtained by [2, Theorems 2.10 and 2.16] in 
the saturated region cannot be achieved. 



52 



We now examine the location of the zeros of the above orthogonal polynomials. Denote 
the zeros of the Airy function Ai by 



> azi > az2 > • • • . 

It is well known [35] that the zeros of Ai lie in (—00, 0). Denote the zeros of pniv, 1/-^) in 
decreasing order yi^N > 1/2, n > ■ ■ ■ > yN,N- We have, 

Theorem 5.7. For fixed k, 

yk,N = b + 2a+ -aik + 0(max + 3 , \b\/N^+-s, l/N's)), (5.45) 

{N^{b + 2a,l/N))l 

and if 6 - 2a < 

yN-k,N = 6 - 2a — -atk + 0(max (l/iV* + 3, |6|/iV^ + ^ 1/iVt)), 

(iV^(6-2a,l/iV))§ 

(5.46) 

Here = q-\{^). 

Proof. From the relation between the equilibrium measure and pi given in Theorem 3.12. 
We see that for t = 1 and fixed N, pi is an increasing function of y and therefore has an 
inverse function. Equations (5.37) or (5.42) and the arguments in [35, p. 406-407] show 
that 

yk,N = b + 2a+ — ^azfc + 0(max \b\/N^+K l/N's)). 

Equation (5.45) can now be obtained using (2.10) and (2.21). Likewise equation (5.46) is 
a consequence of (5.39) and (5.40). □ 
For 6 — 2a > note that is strictly increasing implies that for fixed N,Ti/N{y) 
given by equation (5.29) is strictly increasing and therefore has an inverse. Using the 
notation above we have 
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Theorem 5.8. Suppose that 6 — 2a > then for fixed k, 



1 

yk,N = b + 2a+ aik + 0(max (l/N^ , \b\/N^ , l/N))/N^/\ (5.47) 

{N^{b + 2a,l/N))i 

If I C (0, 6 — 2a) is a closed interval. Then the zeros of pniy, ^/N) in I are given by 

!' = rr;«(^ + <'(i/^))- 



Proof. Equation (5.47) follows from (5.42) as above. To prove the second part note that 
in /, Ai and Bi are nonzero. Furthermore in equation (5.43) the term multiplying the sin 
is exponentially decreasing while the term multiplying the cos is exponentially increasing. 
So that from (5.43), a zero of PN{y, 1/-^), is given by cos{NT ]si {y)) = —Rsiy^N), where 
\Rs{y, N)\ < max{l/N^ ,\b\/N^ ,1/N). □ 

VI. Model Coefficients 

In order apply the above results to the special functions we are interested in we 
now consider in more detail special cases of recurrence coefficients associated with q{t) = 
(t + si) « , Q! > 0, si > so that 

jt/e + si)^ (t/e + si + 1/2)^ 

a(t,e) = a r and b(t,e) = b r- (6.1) 

(l/e + si + l/2)* (l/e + s, + l/2)i 

To make the connection with classical orthogonal polynomials the recurrence coefficients 
given by equation (3.24) can be further divided into the following subcases, 

case 1 6 + 2a > 0, 6 — 2a < (Asymmetric Freud, Miexner-PoUaczek case), 

case la 6 = (symmetric Freud, Hermite case), 

case 2 2a + 6 > 0, 6 — 2a = (Laguerre case), 

case 3 6 + 2a > 0, 6 — 2a > (Discrete Freud, Meixner case) 
The zero distribution for the above cases has been discussed by Kuijlaars and Van Assche 
(see [25] and [26]). 

In what follows we will always take = goinz -^j^ere we use the branch of the log 
that is analytic in a neighborhood of the positive real axis and positive for large positive 
z. From Theorem 3.1 we have the result, 
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Lemma 6.1. With q{t) = (t + si) « and e sufRciently small we Gnd, 

^^^''^-\Ar6(^(i + i)if)"-6fc(«i,-y,6)(l + rr(y,6)) -oo<x<0 
1 1 V 

Q{y,e) = A2e{q{- + -)^)'^-ek{si,y,e){l + r2{y,e)) < x < oo for case 2 



e 2'2a' 
1 1 y 

e) = A3e(Q(- + -) — )"- e/c(si, 2/, e)(l + r2(2/,e)) < x < oo for case 3. 



Here 



with 



k{si,y,e) = - si + - +ln 



= a /\'^-Mn ( - - 1+ W- f- -2 ) Uw. 
Jo \u y w V« . 

For J/ > 0, r^, r2 and have an analytic extensions to C \ (— oo, 0] and obey the bound 

\r+{suy,e)\<Cmax{e^,\b\e^}, |ri(si, y, e)| < Ce^ z = 2, 3 (6.2) 

uniformly on compact subsets of C \ (— oo,0]. In the case of 1 with y < 0, has an 
analytic extension to C \ [0, oo) and obeys the bound 

|rf (6'i,y, e)| < Cmax{e~, \b\e~} 

uniformly on compact subsets of C \ (— oo, 0]. 
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Proof. From the second equation in (3.29) we find for case 1 and y > 0, 

-'"'{1^2) I 



b 

- — + 



2aw 2a 



2aw 2a 



Since y > extracting In(^) from the second term in the above equation yields k{si,y, e) 
and the remaining term can be expanded in a Taylor series to yield the bound (6.2). The 
remaining cases can be derived in a similar manner. □ 

Remark. When 6=0, then A+ = A+ = ^ffj^ . 

We also find that Stirling's formula can be applied to K{n) in Lemma 5.1 Stirlings 
formula to obtain, 

«(n) = ,\ (1 + 0{e)). 

^ ' (27r)2^(si + i)(^i+^)/«^ ^ " 

From formulas for p given by equations (2.2), (2.16) and (2.17) simplify. We find for case 1, 



iV2/3pi(l, l/N) = (AT + + l/2)2/3pi(y), 



(6.3) 



where. 



and 



with. 



Pi(y) = ' 



1^ 

\ V 2a ^ ^ ^TT^ y/4a'2w'2-(y-bw) 



{Ipiyfi^'" I -cosh-y^) 6 + 2a>y 

2/3 



2/3 



(6 + 2a)gi/iv(0) <y<b + 2a, 



iV2/3p2(l, 2/, 1/iV) = (AT + + l/2)2/3p2(y), 



P2(2/) = < 



6-2a 



{bw — y)'^ — Aa^ w'^ 



X 2/3 

cosh-^(^)j y<6-2a 

2/3 



b-2a ^J'ia'^w'^ — {hw—y)'^ 

b-2a<y< {b - 2a)qi/N{0). 



(6.4) 



(6.5) 



(6.6) 
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For case 2, 

Ar2/V(l,y, l/N) = (iV + + l/2)2/3p,(^), (6.7) 



where, 



2/3 



1)1 y 7" - cosh-Vf - 1) ) 4a < y 

-(§)* (--'(1 - - y/i ;«£fe) < ^ S 4a. 

For case 3 and N large enough so that (6 + 2a)gi/jv(0) < 6 — 2a we have 

iV2/3pi(l, l/N) = (iV + Si + l/2)2/3pi(2/), (6.9) 

where, 

(i)M?/ fi'^'^ / -cosh-^^) b + 2a<y 

f -1 

-(1)1 cos-i(^)-2/r\ , ^° b-2a<y<b + 2a 



and 



with. 



(6.10) 

Ar2/3p2(i,2;, i/AT) = (AT + + l/2)2/3p2(y), (6.11) 



-'i'* (^/k V./;-(t-.,. -«'--'(^)) 6-2a<j,<6 + 2<. 

P2(y) = < 



(|)§(cosh-\^)-y/i 



2/3 



Y ^20^* b-2a -y/ — j/)2— 4a2w2 

(6- 2a)Qi/Ar(0) <y<b-2a. 

(6.12) 

For case 3 the contraining measure has a particularly simple density [26] (see equa- 
tion (3.35)) given by 

da q{N + l/2r «_i 



dy N 

where 



Ca2/""\ 2/> (6 + 2a)Qi (0), (6.13) 



a / ^ "'^ (6.14) 

7r76-2a V4a2 - (i^ - 6)2 ^ ' 
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Also 

NT^/N{y) = (iV + si + l/2)7rc«y" - ttsi (6.15) 

Remark. The above formulas can be extended into the complex plane with the use of 
Theorem 3.12 which also gives the relation between the equilibrium measure and p. 
The above results allow us to recast Theorems 5.5 and 5.6. 

Theorem 6.2. Suppose a{t,e) and b{t,e) are given by equation (6.1) with b — 2a < . 
Let ai(t, e), 6i(t, e) satisfy equations (2.24) and (4.20) and set, 

. r(.i + i)ia-^+v^ 

K = J 7-. (6.16) 

(27r)^(iV + si + 1/2)27J 

For y eC+\[b-2a,b + 2a] 

{{N + si + l/2)^yy'pN{y, 1/N) = 

kK,{, ^_^)^e'^(^+^l+V2)/>«-Mn(.-,K.)+V(^-l(-.^)^-l)rf«(l + <^^(^,l/^^)). 

[y — by — 4a^ 

(6.17) 

where zi («, y) = ^ . The bound 



|(/>i(y, 1/N)\ < max(6/iVV«, l/iV^/«, l/N), (6.18) 

holds on compact subsets of C+ \ [b — 2a, b + 2a\. For y e (0) given by Lemma 3.10 and 
N sufRciently large, 

{{N + si + l/2)iyy^PN{y,l/N) ^ 

g(JV+si + l/2)A(^)« 

20r/««i / pi{y) ^ 



(iV + si + l/2)2Tr-i \{y-by-4a'^ 



{Ai{{N + si + l/2)2/3pi(y)) + ri(y, l/AT)) 

(6.19) 

where pi is tiie extension of (6.4) for 6 — 2a < or (6.8) if b — 2a = and A = A'^ for 

b — 2a < or A2 for b — 2a — 0. The error term ri satisfies, 



|ri(2/,l/iV)e^(^+^i+^/2)^?(^'^/^)| <dmax(|6|e",e«,e), (6.20) 
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uniformly on compact subsets of Lf{0). Ifb — 2a<0 and y e -£'2^(0) given by Lemma 3.10 
then for N sufRciently large, 

(_l)W((Ar + s, + l/2)ii-y)y^pN{y, 1/N) _ 

g(iV+si + l/2)A-(-^)« 

(6.21) 

wiere p2 tie extension of (6.6). Tie error term r2 satisfies the bound, 

|r2(y,l/iV)ei(^+^i+^/2)p|(j/,i/Ar)| < rfmax(|6|e« , e« , e), (6.22) 



uniformly on compact subsets of -£'^(0). 

Likewise for 6 — 2a > we find, 

Theorem 6.3. Suppose a(t, e) and b{t,e) are given by equation (6.1) witi b — 2a > . 
Let ai(t, e), bi{t, e) satisfy equations (2.24) and (4.20). For y e C+\[0,b + 2a], 

{{N + si + l/2)iy)>Ar(y, 1/N) = 

(6.23) 

wiere 5i(u, j/) = 2^:^ ~ ^ ^^<^ ^■s giVeii by (6.16). Tie bound, 

101 (y, l/iV)| < max(6/iVi/«, l/iV2/«, 1/iV), (6.24) 

holds on compact subsets of C+ \ [0, 6 + 2a]. For y e -^5*^(0) giyen bj Lemma 3.11 and A?" 
sufficiently large. 



{{N + si + l/2)^yY-pN{yMN) ^ 

\ 1/4 

(Ai((Ar + si + l/2)2/3pi(y)) + ri(l, y, l/N)). 



g(Ar+si+i/2)(^)« 

2^Rki ( pi{y) ^^^"^ 



(iV + si + 1/2)2^-^ V(^-^)^-4a2 

(6.25) 
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where p\ is the extension of (6.10). The error term ri satisfies (6.24) uniformly on compact 
subsets ofL^iO). For y e {0,b + 2a) 

{{N + s, + l/2)iyy^pr,{y, 1/N) ^ 2^kK^ ( \ ^'^ 



(6.26) 



(6.27) 



g(iv+.i+i/2)(^)« V > (iV + si + l/2)2fe-H4a2-(y-6)2 

(sin7r((iV + si + l/2)y«c« - si)r + i?i(y, 1/iV)) 

X (Ai((iV + + l/2)2/3p2(2;)) + ri(y, 1/iV)) 
+ (cos7r((iV + si + l/2)y"ca - si) + i?2(2/, 1/iV)) 

X (Bi((iV + 51 + l/2)2/3p2(2/)) + r2(|/, 1/iV)) 

where uniformly on compact subintervals of (0, 6 + 2a) 

|fi(|/, i/iv)ei(^+*i+^/^)^2^(^'^/^)|, |f2(2/, i/iV)e-|(^+*i+V2)4(J/.V^)| 
|i?i(y, l/iV)|, ||i?2(y, 1/A^)| < rfmax(|6|/Ari, l/AT^, l/AT). 

Remark. In the above cases ^(6 + 2a, 1/N) = °'^^'\^^^2a)'^^ so that the formulas for the 
zero's given in Theorems 5.7 and 5.8 simphfy. 

VII. Applications 

Example 1. Hermite polynomials, 6 = 

We consider well known classes of orthonormal polynomials associated with the various 
cases given in the previous section. We begin with the orthonormal Hermite polynomials 
{Hri{x)} which have 

w{x) = e~^^ 

as their weight function (35, 41). In this case bn = 0, n > 0, an = a/¥, n > 1 and 



Qn = SO that a(t, 1/A^) = ^ n+i'/2 - Theorems 3.13 and 4.3 we take ai[n) = a[n), 
bi{n) = so that ki = 1 and Ai = 1. Since po{y) = 1, Hn = -174- In order to compare 
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the previous results with the hterature we make the scahng y — > -\/2j/, p\{y) = — ^-^Pi{y) 
and Ajv = -\/2A/" + 1 so that x = XnV- Thus 



(l (yVy"^ - 1 - cosh ^ y)) ^ for 1 < y 
(^f (^cos-iy-yVl -2/^))^^^ for (0) < y < 1. 



Pi{y) = 

[ (f(cos 

Since si — 0, b — 0, a — 2 and = 1 we find uniformly on compact subsets of C \ [—1, 1] 



^^(^^^) = i^^5!z!r!e-*^v/^ (1 + 0(1/N)) (7 1) 



Also on L^(0) 

1/4 



-f^jv(Ajvy) _ \/2 / pi(?/) 

1/ 

AT 

where uniformly on compact subsets of -£'i"(0), 



3 

|el^'^p?(s')ri(y, 1/7V)| < d/N. 
On L+(0) 

where uniformly on compact subsets of L^(0), 

3 

|e§^^p?(-f)r2(2/, 1/^)1 < d/N. 



Formula (7.1) was first derived from the point of view of difference equations by [42] and 
formula (7.2) by [14]. For y > this formula is found in Olver [35, p. 493] (see also Szego 
[41, Thm 8.22b]). From Theorem 5.7 we find that the zeros of i/'jv(Ajvy) are given by, 

yk,N = 1 + + 0{1/N^), yN-k,N = -yk,N- 



(See also Szego [41, p. 132].) 



61 



Example 2. Meixner-PoUaczek polynomials, 6 + 2a > 0, 6 — 2a < 

The Meixner-PoUaczek (Meixner polynomials of the second kind) Mn{x,d,r]) are or- 
thonormal with respect to the weight function 

Since po{x) = 1 we find that [22] 



The coefficients in the recurrence formula for the orthonormal Meixner-PoUaczek polyno- 
mials [6], Mn{x '.6,7]) are given by, 

ai(n) = V(^2 + lV:^(n + r/-l),r/ > 0,5 > 0, 

and 

6i(n)= (n+l) 2S 
In this case we will use a comparison system whose coefficients are given by, 

a(n) = ^/ S"^ + 1 -|- ' ^l*^) — ^i(^)) 

so that 



With q{t) = t + ^ we find si = a(t, e) = a ^^^/|\^ , e) = 6 ^|f+!j = e) and 
ai(i, e) = ^'llVh r^ where a = VS'^ + 1 and b = 2S. For t strictly greater than zero 
|ai(t, e) — a(t, e)| = O(e^) and sup(g[o,t/i] ""^^(^^ ^)l — 0{e^)- Thus equations (2.24) 

and (4.20) are satisfied. With the scalings y — > 2y, p = (^)^p, Xn = 2N + r] and x — XnV 
we find, 

du , / y(id + 1) 

I in 



\/y'^ - 25yu - \ {Sy + l)i + ^Jy^ - 25y — 1 

= 2tan"^ 



(7.4) 

1 



.y 



+ -25y-l 
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so that, 

/ ln(5i {y, u) + ^/zi{y,uY - l)d 
Jo 



= cosh ^ (- — -]+iy In' + ) 



Also 



{Sy + l)i+y/y^-26y-l 



-(!C0S-1(^) - f|/COsh-l(f + A))2/3 fo, ^^/^(O) < < 1 



which can be extended to -£'^"(0). We also find that A'^ = a cot ^{5). In a similar manner 
— a(7r — cot~^(5)) and 

(l(-y)cos-i(-f - 4) - I cosh-i(^))2/3 for < 1 

-(|C0S-H^) - t(-y)cOSh-^(-f - i))V3 fo, ^^/^(o) < < 1, 

which can be extended to Lj(0). Thus uniformly on compact subsets of C\[S — a,S + a], 



P2{y) = 



+ V2/'-25y-l)^ 



7r{XN)Hy^ - 25y - 1)4 \{5y ^ l)i + -25y-A) 



(7.5) 



On L+(0) 

— 1 ^ 

{\Ny)~MN{\Ny ■■ 5, ??) 



gCot-i(<5)Ajvy/2 

' P'^y' ^ {Ai{Xyp,{y)) + r,{y,l/N)), 



(7.6) 



^X]i' W-25y-l 
where uniformly on compact subsets of L^(0), 



3 

e§A^p?(3/)ri(y,l/iV)| < d/N. 



On i^J(O) we have, 



g-(7r-cot-i(5))AArJ//2 

/ ~ / X N 1/4 (7-7) 
' P'^y^ ^ (Ai(A^/V~2(2/))+r2(2/,l/iV)), 
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where uniformly on compact subsets of -£'2^(0), 

3 

|e§^'^p2^(j/)r2(y, 1/N)\ < d/N. 



From Theorem 5.7 we have for k fixed, 



yk,N = S+VP^+l ^^f^ -] (^^^1 azfc + 0(l/iVS), (7.8) 



A 



TV 



and 



TO-t,N = *- V*" + l- I I 1^-^; I ait + 0{l/m). (7.9) 

Equation (7.5) was first obtained by Geronimo and Van Assche [18] with a slightly different 
scaling. With the use of equation (7.4) this is similar to the asymptotics obtained by Li 
and Wong [30] equations (6.9) and (6.10) with the scaling X^y = 2Ny. Equations (7.6) 
and (7.7) agree with equations (6.23) and (6.24) of [30] when the above scaling is taken into 
account. The location of the zeros given by equations (7.8) and (7.9) should be compared 
with equations (7.17) and (7.18) of [30] (see also Chen and Ismail [5]). 

Example 3. Laguerre polynomials, 6 — 2a = 

The orthonormal Laguerre polynomials L^{x) have weight function 

w{x)—x"e~'^, x>0, a>—l, 

with recurrence coefficients, 

ai(n) = \/ n{n + a), and bi{n) = 2 H ^- — ^ 

Thus L^(x) = _£si£)=. As a comparison system we will use, 

a{n) = n + ^, h{n) = hn- 



With q{t) = t + f we find si = f , a(t, e) = ^^i^, 6(t, e) = = ^>i(^,e) and 

/ — ( — -\-oc) 

ai(t, e) = -/iliN where 6=2. For t strictly greater than zero |ai(t, e) — a(t, e)| = 

64 



O(e^) and sup^g[o,t/i] " = 0{e^). Thus equations (2.24) and (4.20) are 

satisfied. To move the oscillatory interval to [0, 1] we scale y — > 4y, and set p = 4rsp, 
Xn — 4:N + 2a + 2 and x — XnV- In this case A2 — 1, 

_ yr(^Ti) 
r(f + i) ' 

and 

^'^^^ l-(|(|cos-H2y-l)-^/^^))^/3 for 0<y<l, 

which can be extended to L^{0). Thus uniformly on compact subsets of C \ [0, 1] 



(Aiv2/)^^«_ 1 2y-l + 2^ 



4 



V2^iy^-y)-^\ e^x/^ 

On i^]'"(0) we have 



{1 + 0{1/N)). (7.10) 



'^"^''i;'^''^' = ^ (^) ' (^(4hiy)) + i/iV)). (7.11) 



where uniformly on compact subsets of -£'i"(0), 



3 



|e3^^Pi (s')ri(y,l/Ar)| < d/AT. 

Starting from the recurrence formula equation (7.10) with a different scaling was first 
given in [18]. Equation (7.10) can be found in Szego [41, Thm 8.22.8b] for y > 1. See also 
Frenzen and Wong [12]. We have for the zeros of L^(Ajvy), 

4I 

yk,N = 1 + ^aik + 0{1/nI). 

(See [41, p. 132].) 

Example 4. Meixner polynomials, 6 — 2a > 
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For the orthonormal Meixner polynomials rhn{x : (3,c) have an atomic measure with 
masses 

w[x) = ^c', x = 0,l,..., 



and recurrence coeflBcients 

1 -I- /- / r \ 

<c<l, 0</?. 

1 — c 1 — c\ i-|-cy 

Thus 



ai(n) = -^Vn(n + /?-!), and 6i(n) = + y^/^I ' 



mn{x : p, c) - 



(l-c)f 

We apply the scaling x — > ^ej/ — Ai, where Ai = | and is given below. The shifted 
Meixner polynomials rhn {x : P, c) satisfy the recurrence formula 

ai{n+ l)m„+i(5) + {bi(n) - x)-fhn{x) + ai{n)rhn-i{x) = 

where ai{n) are as above and hi{n) = jz^in + f )• 
The comparison system is 

a{n) ~ — ) ' ^^^^ ~ bi{n). 

Thus 

_ r(/3)V^ 

With q{t) =t + si,si^ ^ ait, e) - a^^i^, e) = ^||tjfy = e) and ai{t, e) = 
a ^^||^f -y^ where a = and 6 = we see that ^ > 1- For t strictly greater than 
zero \ai{t,e) — a{t, e)| = O(e^) and sup^^jQ^^^.j |af (t, e) — a^(t, e)| = O(e^) so equations (2.24) 
and (4.20) are satisfied. In this case 



pi{y) - 

and 

h{v) 



(f)l(ycosh-i(A __!_)_ cosh-^(i^))2/3 for 6 + 2a < 

2 



-(§ )i (cos-i(^) - y cos-i(A - ^))2/3 for 6 - 2a < y < 6 + 2a, 



-d )§(y cos-i(^ - A) _ eos-i(^))2/3 for 6 - 2a < y < 6 + 2a 
(f )i(cosh-^(^) - ycosh-^(^ - A))2/3 for < y < 6 - 2a. 
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We also find that = 2aln(^) and from (6.14) ci = 1. With Ajv = + /?/2, we find 
uniformly on compact subsets of C \ [0, 6 + 2a] 



(Aiv2/) 2 - 



(ttl\XNy 
y 2a ) 

r(/9) ^c^ {{y-h+ y/y^ - 2by + l)/2a)^^ 
v/2^(y2-26y+l)3 



2ay 



6y - 1 + - 26y + 1 



(l + 0(l/iV)). 
(7.12) 



On [h — 2a, oo), 



V 2a I 



2a J - (Aiv)* Vy'-2&y+l 

where uniformly on closed intervals of [h — 2a, oo). 



(Ai(A^pi(y))+ri(y,l/iV)), 



3 

ei^^^?(«Vi(y,l/iV)| < d/N. 



On (0,6 + 2a), 



(Aiv|/)^mjv(Ajvy- § :/3,c) ^ ._..iv\/2r(^ 



1 £ 

2C4 



(k±l\\Ny 
\ 2a I 



Pliy) 



COS 



y2 _ 26y + 1 

s f TT (Xj,y -^)+ R,iy, l/N)) (Ai((A^/'p2(2/)) + ^(z/, 1/A^)) 



+ ( siuTT ( A^y - ^ ) + R2{y, l/N) ) (Bi(A^/=^p2(2/)) + ^2(2/, 1/iV)) 



where uniformly on closed intervals of (0,6 + 2a), 

\ei^^p|^y\^{y,l/N)l le-i^-^^^^^Vsd/, l/iV)|, \R,{y,l/N)l \R2{y,l/N)\ < d/N. 



From Theorem 5.8 we find 



yk,N 



1 r 



+ 



1 -c 



2 

A^ 



azfc + 0(l/iV3). 



N 



and 



yN-k,k 



{2k-l)7r + p 



2A 



JV 



Example 5. Continuous dual Hahn polynomials 6 — 2a = 
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In order to keep the exposition contained we will assume the parameters ai,6i, and 
ci are positive and ai + 61 + ci > 3/2. In this case the normalized continuous daul Hahn 
polynomials hn{x; ai, bi, ci) are orthonormal with respect to the weight [22], 



wix) = i- 



r(ai + ix)T{bi + ix)T{ci + ix) 



2 



a; > 0. 



T{2ix) 
If 

An ^ {n + ai+ 61) (n + ai + ci) and C„ = n(n + 61 + ci - 1), 
then the recurrence coeflBcients are 

bi{n) = An + Cn - a-i^ = 2 [ n + (ai + 61 + ci - l/2)n H 



and 

al{n) = An-iCn = n{n + ai + 61 - l)(n + ai + ci - l)(n + ci + 61 - 1). 
In this case 

hn{x;ai,bi,ci) = 



v/r(ai + 6i)r(ai + ci)r(6i + ci)" 
We chose the comparison system with recurrence coefficients, 

a{n) = q{n) and b{n) = 2q{n+ 1/2), 

where 

q(n) = {n+ (ai + 61 + ci - 3/2)/2f. 

Thus 

^ v/r(ai + bi)r(ai + ci)r(6i 



4 



With m as above = 2ai+2b,^+2ci-3 ^ ^ ^ ^ _ai(j|_^ 

b{t,e) = 2-^^^T^ and bi{t,e) = J'(]!'^l^ we find for t strictly greater than zero \ai{t,e) — 
a{t,e)\ = 0{e^) and sup^g[o,t^.] |af(t,e) - a'^{t,e)\ = 0{e'^) = sup^g^.t^i] \h{t,e) - b{t,e)\ 
so equations (2.24) and (4.20) are satisfied. We make the scaling y ^ Ay and set Ajv = 
{2N + ai + 61 + ci - 1/2)2. For y e C \ [0, 1], we find that 



/' 

^0 



^Jwy-w^ \ \y J J ^ \ y y 
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so that 



/ tt 2 In 
Jo 



(7.13) 



Also 



2 

2 



^ ^ (^cos-^(^ - 1) - cosh"^(2|/ - for 1 < y 



-p(cos-i(l - 2y) - Vycosh-^(| - for g^(0) < y < 1, 

which can be extended to L^{0). Thus uniformly on compact subsets of C \ [0, 1], 

/ 2oi+26i+2ci — 3 s 

{XNVr ^ 'hniXNy;ai,bi,ci) = 
(2y - 1 + 2 ) ^ (l - 1 + 



V27rA^(y(y-l))^ 
For y e L+(0) we have, 

(Xnu)^ 4"^ ' ^/in(Ajvy;ai,6i,ci) _ 
— ]_ 

^^(j|^)'(A'(4p.fe)) + n(y.i/iV)). 

where pi = -^jzPx- Uniformly on compact subsets of -£'^"(0) 

|e§^^^?(2')ri(y,l/Ar)| < d/iV. 



-(1 + 0(1/Ar)). 



From Theorem 5.7 we find 



y.,iv = i + ^ + 0(i/ivt). 



Example 6. Wilson polynomials 6 — 2a = 
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The Wilson polynomials are the most general hypergeometric polynomials in the 
Askey scheme which do not have a purely discreet orthogonality measure. We denote 
the orthonormal Wilson polynomials by Wn{x'^ '■ ai, 61, ci, (ii) and for simplicity wc will 
assume that Re(ai, 61, ci, (ii) > 0, non-real parameters occur in conjugate pairs and 
ai + bi + ci + di > 1. These polynomials are orthonormal [22] with respect to the weight 

2 



w{x) — 



r(ai + ix)T{bi + ix)T{ci + ix)T(di + ix) 



V{2ix) 



, X > 0. 



With 



and 



^ _ (n + ai + 61 + ci + (ii - l)(n + ai + hi){n + ai + ci)(n + ai + di) 
(2n + ai + 61 + C + 1 + c/i - l)(2n + ai + 61 + ci + 

n(n + 61 + ci — l)(n + 61 + (ii — l)(n + ci + (ii — 1) 



Cn 



(2n + ai + 61 + ci + di - l)(2n + ai + 61 + ci + di - 2) 



we have 



In this case, 



ai(n) = •\/^n-iC'n, and fei(n) = ^„ + Cn - . 



2 7 7 N / r(ai + 61 + ci + rfi) 

W^(a; ■aiM,cudi) = \ — — — - — r — ——Pn[x). 

r(ai + Oi)r(ai + ci) ■ ■■V{ci + di) 

We chose the comparison system to have recurrence coefficients 

= and h{n) = 1/2), 

where 

, , / ai + &i + ci + cZi - 2^ 
q{n) = ( n + ^ 

This comparison system is very similar to the one used for the continuous dual Hahn 
polynomials so those results can be used here. With q{t) as above si = ai+^i+ci+c'i-^ ^ 
A, = ^, a{t,e) = i^, ax(M) = KM) = and 6i(M) = ^ we 



find for t strictly greater than zero |ai(t, e) — a(t, e)| = 0{e ) and sup^g[o,tji] kil^?^) ~ 
a^(t, e)| = O(e^) = sup^g[o,t/i] l^i(^) ^) ~ ^(^) ^)l so equations (2.24) and (4.20) are satisfied. 



70 



We set Aat = (iV + «i+k+ci+di^)2_ ^sing equation (7.13) we find that uniformly on 
compact subsets of C \ [0, 1], 

{XNVr ^ 'WniiXNy) ; ai,bi,ci, di) = 

{2y-l + 2.M^)fT{l - I + ^2Jl^))^ 
, ^ ^— ^ (1 + 0(l/iV)), 

while for y e -£'i'(0) equation (7.14) gives, 

(Ajvy)^"'^"'^^'^"' '^W^n((Ajvy)^; ai, bi, ci, di) ^ 

gTT-y/AjvJ/ 

Pl(?/) 



V2^A* Vy(y-i) 



(Ai(A|,pi(2/)) + ri(2/,l/iV)), 



AT 



where uniformly on compact subsets of L^(0), 

\el^Pi^y)ri{y,l/N) \ < d/N. 

From Theorem 5.7 we find 

yk,N = l + ^ + 0{l/N-s). 

Appendix A. 

In this appendix we sketch a proof of Theorem 3.13. We begin by considering solutions 
to the differential equation 

-^ = e px. (A.l) 

We assume that p is an infinitely differentiable function of t > The above differential 
equation shows that the Taylor expansion for x(e~3p(t ± e)) can be written interms of 

2 2 

x(e~3p(/;)) and x(e~3p(t)). We examine this expansion by setting, 

where Rn is the remainder term in the Taylor expansion. Wang and Wong [44] essentially 
show, 
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Lemma Al. Suppose x is any solution of (A.l) then, 

= xi^~hit))Xi{p, p±, ±e) ± e^x'(e~^PW)^2(p, P±, ±e), 

and 



X'{e sp(t)±eBp±) = 

± e-3x(e-ip(t))X3(p, p±, ±e) + x'(e"^p(t))X4(p, p±, ±e), 



where 



and 



If 



then for p > 0, 



and 



1 F± 

Xi,k{p,P±) = ±^ / sXi,fc-i(p,s)sinh^(p± - s)ds 

VP Jo 

/ sXi^k-i{p, sp±) smb. ^/pp±{l - s)ds. 
Jo 



(A.2) 



(A.3) 



X3 = ^Ai, (A4) 



X4 = ^X2. (A.5). 

ap 



X,(p,p±,±e) = J](±e)-X,,,(p,p±) z = 1,2, (A.6) 

n=0 



smh. l~oo 

^1,0 (p,p) = cosh^p, X2,o (p,p) = — — — , (A.7) 



(A.8) 



Above if p < then is replaced by iy/—p- 

If (A.2) is differentiated twice with respect to p±, then the Airy differential equation, 
and the conditions 

Xi(p,0,±e) = 1, ^ = 0, 

dp± 



and 



dp± 
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yield upon equating powers of e the following differential equations [44] , 

- pX,,o, Xi,o(p, 0) - 1, - 0' 



op± dp 



± 



and 



dpi ^^^'''^ X2,o(p,0) = 0, ^^^ = 1' 

= pA2,fc ± P±A2,fc-l A2,fc(p, U) = U, ^ 

OP± dp± 

Integrating these equations yields (A. 7) and (A. 8). 

Using induction Wang and Wong essentially show that, 

|3fc 



and 

Substituting the formula for p± into Xj^„(p, p±) yields 

m 

^i,n(P,P±) = ^{±(^y Xi^njiP, P) + Rm, 

Following [44] and [14] we now look for solutions of the equation 

a{t + e)i;{t + e) + {b{t, e) - y)i;{t, e) + a{t, e)i;{t - e) = 0. 



Write, 



n=0 n=0 



and 



n=2 
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The assumption on the form of b{t) allows us to seek solutions of the form 

m = x{e-h{t))C{t, e) + e-^x'{e-h{t))D{t, e), (A.16) 

where 



oo 



C{t,e) ^J2^s{t)e' and D{t,e) = J]L>,(t)e*. (A.17) 

Substitution of equations (A.16), (A. 2) and (A. 3) into equation (A. 12) then equating 
coefficients of x ^^d x' yields 

a{t + e)(Xi(p, p+, e)C{t + e) + eX^ip, p+, e)D{t + e)) + {b{t) - y)C{t) 

(A.18) 

a(t)(Xi(p, p_, -e)C{t - e) - eX^{p, p_, -e)D(t - e)), 

and 

a{t + e)(X2(p, p+, e)C(t + e) + eX4(p, p+, e)D{t + e)) + e(6(t) - y)D{t) 

(A.19) 

a(t)(-X2(p, p_, -e)C{t - e) + eX4(p, p_, -e)L>(t - e)). 

Substitute the expansions for C, D, o, 6 into equations (A.18) and (A.19) then collect 
powers of e and use equation (A. 8). For e'^ we find from (A.18) 

(2aoXi,o,o + 6o - y)Co = 0, (A.20) 

Thus (A. 7) shows 

Xi,o,o = cosh( Vpp') = (A.21) 

zao 

For the coefficient of e in (A.19) we find 

2ao[(X2,o,i + X2,i,o)Co + X2,o,oC'q] + aiX2,o,oCo + (2aoX4,o,o + {ho - y))Do = 0. 

Equations (A. 5) and (A.20) show that the term multiplying Dq is equal to zero and (A. 7) 

gives, 

sinh ./pp' p" _ , , , 

^2,0,0 = , ^2,1,0 = ^ cosh Vpp', (A.22) 

and 

^2,0,1 = ^ cosh Vpp' - ^ sinh Vpp'. (A.23) 
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Thus 



Co 2ao 2sinh^p'l2Vp "^^^ i 4p' 



which has Co{t) = 5r(t) as a solution. To continue on we utihze (A. 4) to obtain 



^3,0,0 = Vpsinh^p', X3,o,i = ^cosh ^pp', 



and 



-^^3,1,0 = 1^ sinh y/pp' + cosh ^/pp' . 



If we now look at the coefficient of and use the results just obtained we find by induction 
that, 

2ao(X3^o,i + -'^3,i,o)-Dn-2+ai-'^3,o,o-Dn-2 + 2aoX3,o,0-D^-2 

— 9n{Co, . . . , Cn-2, -Do, • • • , Dn-4:), n > 2, 



which can be recast as 



d r 
di 



{p4alsmh\^p')y^^ 

1/4 

gn{Co, C2, . . . , Cn-2, Dq, . . . , Dn-4:), n>2. 



1 



Likewise 

d 
di 



p4a§sinh {^p') 



al sinh^(ypp') 



1/4 



a, 



n-2 



3/4 



agsinh {^p') 



/n(Co, C2, . . . , Cn-3, Dq, . . . , D^-s), n>2 



(A.24) 



(A.25) 



where fn and Qn are linear in the variables in their arguments and their higher derivatives. 
In the above equations we take the convention that Dn = = if n is negative. 
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Theorem A2. Suppose (3.24) and ic)-iuc) hold, bo = b{t,e) and bi{t,e) = 0,i > in 
equation (A.15). In equations (3.52)-(3.54) let p ^ = {t+ - t)(Ci)^/^ and let Li(e) be 
given as in Lemmas 3.10 or 3.11 and L^{e) = Li{e) fl C-^-. If y E L^{0) then there exists 
an ei such that y e -^i^(e), < e < ei. Let x be any entire function solution of (A.l) and 
ip = gx- For each e, < e < ei and all t : t e [tin, tfi] 

a{t + e, e)'ip{t + e, y, e) + {b{t, e) - y)ip{t, y, e) + a{t, e)ip{t - e, y, e) = y, e), (A.26) 

where 

\l3{t,y,e)\<c{y,e)e^ sup [|x(e"3p(«, y, e))| + |x'(e"^p(«, y, e))|] . (A.27) 

If K is a compact set in L^{0) then there exists an such that c{y, e) is uniformly bounded 
onKx [0, ex]. Furthermore for fixed {t, e), /3 e H{K nC+) and f3{t, •, •) e C{K x (0, e^]). 
Ifb-2a<01etp = p2 = {t~ - t){C2f/'^ or ifb-2a> let p = p^ = {t - t-){C2f''^ in 
equations (3.52)-(3.54). For both cases let L2{e) be given as in Lemmas 3.10 or 3.11 and 
Lj(e) = L2{e)r\C+ .Ifye -£'^(0) then there exists an ei such that y e -£'2^(e), < e < ei. 
Let X be any entire function solution of (A.l) and ip = gx- For each e, < e < ei and all 

t . t & [tin, tfi] 

a{t + e, e)'ip{t + e, y, e) + (y - b{t, €))'ip{t, y, e) + a{t, €)'ip{t - e, y, e) = /9(t, y, e), 

wiere the error (3 satisfies (A.27). If K is a compact set in L'^{Qi) then there exists an 
ex such that c{y,e) is uniformly bounded on K x [0,6^]- Furthermore for fixed {t,e), 
(3 e H{K n C+) and P{t, ■, •) G C{K x (0, e^])- 



Proof. We first take y real and choose ei < eo so that tin — ei > 0. From equations (A.7)- 
(A.ll), (A.21)-(A.23), Theorem 3.12 and the properties of q we find 

|Xi,o(p, P±, ±e) - (Xi,o,o ± e^i,o,i)| < c(y)e^ z = 1, 2 (A.28) 
|X,,i(p, p±, ±e) - Xi,i,o| < c(y)e, z = 1, 2 (A.29) 
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and 



where c{y) is independent of (t, e) e [Umtfi] x [0,ei]. This imphes from (A. 2) that 

1 i-l 



1 1-1 

± eh'{e-ip{t)) E E(±^)'^'^2,^, + e^{t, y, e). 



«=0 7=0 

With the above estimates we find that, 

\eiit,y,e)\ < e'^c{y) sup |x(e"^p(w, y, e))| + |x'(e"^p(w,y, e))| 

ue[t-e,t+e] L 

With the choice of Cq = ^(t), the terms of e*, z = 0, . . . , 1 vanish in equations (A. 18) and 
(A. 19). The remaining terms times x or x' when added together give For y e -£'i'(0) 
differentiation with respect to t of (3.38) and Theorem 3.12 show that cosh(^p') = 
. Thus for fixed (t, e) Morera's Theorem implies the above equations for Xi^k can be 
analytically continued to Lf{0). Theorem 3.12 also shows that for fixed (y,e) in L^{0) x 
[0,ei], ipi and V'2 satisfy equation (A. 26) with P satisfying the bounds and smoothness 
given above. The other cases can be argued in a similar manner. □ 

Remark. Integration of Equations (A. 24) and (A. 25) as in [35, p. 409] allows the extension 
of the above results. However since this is not needed for the applications considered above. 
We now have. 

Theorem A3. Suppose (3.24) and ic)-iiic) hold. Let ai{t, e),bi{t, e) e C([0, oo) x [0,eo]) 
satisfy (2.24) with 6i(t, e) real and ai(t, e) strictly positive on every compact subset of 
(0, oo] x [0, eo]. In equations (3.52)-(3.54) let p = = {t+ - t)(Ci)^/^ and let Li(e) be 
given as in Lemmas 3.10 or 3.11 and L^{e) = Li(e) ("1(7+. Suppose K E Lf{0), K compact 
then there exists an ex such that for each (y, e) E K x (0, e^] and all n : ne E [tin,tfi] 
there exist solutions fi i = 1,2 of equation (2.23) such that 



fi{n) = il)i{n)^-ri{n), 
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(A.31) 



where 

I ri{n) fi{n) - V'i(ri) 



<d{K)e, 1 = 1,2. (A.32) 



Furthermore for fixed e ^Ji^ e H{K n C+) and G C(K x [0, e^]). If 6 - 2a < 

iet p = p2 = (tp - t){C2f/'^ or ifb-2a>0 let p = p2 = {t - t~){C2)^/^ in equations 
(3.52)-(3.54). For both cases let L2{e) be given as in Lemmas 3.10 or 3.11 and I^J(e) = 
L2(e) n C*-!-. Suppose K C LJ(0), K compact, then there exists an sucii that for 
each (y, e) e K x (0, e^] and all n : ne E [tin,tfi\ there exists solutions fi, i = 1,2 of 
equation (2.23) such that 

{-irMn)=Mn)+r,{n), (A.33) 



where 



<d{K)e, i=l,2. (A.34) 



Furthermore for fixed e, e H{K n C+) and e C(K x [0, e^^]). 

Proof. The proof follows closely that of Theorem 4.4 in [14] . Given K a compact subset 
of iv]'~(0) there exists and such that K C for e G [0, e^]- From Theorem A3, 

equations (3.52) and (3.53) we find for j = 1,2 that V'j and u^^^ satisfy equation (A. 26) 
with error j3j and j3^^^ respectively. Also Lemmas 3.10 or 3.11 tells us that for fixed y and e 
Repf {t, y, e) is a decreasing function of t. With tn = en, aj = ^^^^fe^y^, V'j("') = ^§T^i 

^-^A — _ -.T.^ -Rrn 



= P^'^ = and A{w{n)) = w{n) - w^{n), we find using (2.23) that 



«(2)(n - 1) (a^in) - a^^n - l))-u('\n + 1) ""^^^(^,1^)' '^ ^2(n + 1) - a^in)) ^^^^^^ 

= h2{n) + q2{n) 

and 

u^'\n + l){a^{n + 1) - a^{n))-^^^^^^y'\n - 1) (ai(n) - a^{n - 1)) 

a((?i + l)e, e) (A. 36) 

= hi{n) + qi{n). 
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where 



and 



+ (3j{nt 



a((n + 2 - j)e,e)- 



a((n + 2 - j)e,e)' 

Selecting a solution /2 such that (T2{nfi) = ^ = a2{nfi — l) where n/j is the largest integer 
such that Ufie < tfi yields 

where 

a(ze,e)zx(^)(z-l)z.(^)(z) 

The above formula for o can be recast as 

n-fi — l rifi 

a2{n)^ ^ g2(n,0 ^(2y_^^ + E ^2(n, i)a2(i), (A.37) 

i=n+l i=n+l 

where 

K2{n,i) = G2{n,i) A — — + — — - -j—. — — 

+ G2{n, ^ - 1)A ^^^((— T)^ } ^^{[^y 

Here G2{n,i) = G2{n,i) for z < n/j — 1 and zero otherwise. Note that G2{n,i) = for 
i < n. For j = 1 we apply a similar analysis to equation (A. 36) with the conditions 
ci('^m) = = (Ti{nin + 1) where riin is the smallest integer such that > Un yields 

n— 1 / n— 1 

i=ni + l ^ i=ni 
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Here Gi (n, i) = G\ (n, i) for z > ni + 1 and zero otherwise, 



n-l 



a((i + l)e,e)^W(z + l)7i(i)(z) 



k=i 

and 

Ki(n,i) = -Gi(n,i) A — ^ ' \ + — f — ^ ^ , 



+ Gi(n,z + 1)A 

Let 



((z + l)e,e)\ 



a(z + 2)e,e) J u(^'){i + 2)' 



^(t) = |eii^(*'2''^)^|. (A.39) 
We interrupt the proof of Theorem A3 to show, 

Lemma A4. Suppose that (3.24) and ic)- iiic) hold. Let K be a compact subset of -£'i'(0). 
Then there are constants d{K) and so that for all (n, i) : [ne, ie] C [tin, tfi\i^ ^ [0; ^k] 
the following inequality holds, 

\Gi{n,i)\ <Gi{i) =d{K)e-^\u^^\i)u^^\i + l)\, (A.40) 

and 

\G2{n,i)\ <G2(^) ^d{K)e-^\u^^^(^)u^^\i-l)\. (A.41) 
The above inequalities imply that 

J2 Gj{i)<d{K)e-\ J = 1,2. (A.42) 



Proof. For fixed y and e we write (see [14, eqs (3.45) and (3.46)]) as, 
1 

__ 1 r u^^\k + l) _ u^^\k) \ 

" a{k + 1) w(i)(fc + l)u(^){k) - u(^){k)u(^\k + 1) \u(^Hk + 1) w(2)(fc) ) 
_ {k + l)^(^) jk + l)- {k + l)^(^) {k + 1) f u^^\k + l) _ u^^\k) \ 
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where Theorem A2 has been used to obtain the last equahty. Here 



Summation by parts yields 
1 



i-l 

E 

k=n 

i-2 



a(/c+l)w(2)(/c)w(2)(A; + l) 



+ 



«(2)(n); a(z) «(i)(z)«(2) (i-l) -l)«(2)(^)• 

Since -^]'"(0) e S'o U we find from the properties of Aij, i = 0, 1 and w^^\ ([35, p. 238 
and p. 418] and [14]), the nonvanishing of g and the monotonicity of Re pa that for m < i, 



u 



(2)(m) 



< d{K)E-^{if 



Also Theorem A3 implies that 



\(3^'^\k + l)w(^)(/c+ 1)1, \(3^^\k + l)w(2)(/c+ 1)1 < d{K)e^. 

Lemma Al, the Wronskian for Aii and w^^^ and the nonvanishing of g and a{t,e) show 
that 

<d{K)e-^. 



\Du 



In the above inequalities d{K) is independent of e e [0,eK], n : ne e {tmitfi] and y & K. 
Utilizing these inequalities the definition of u^'^\ i = 1,2 and the fact that there are at 
most Ufi e terms in the summation gives equation (A. 40). This bound shows that 

i=nin + l i=nin 
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The above sum can be rewritten as, 

J2 \u^^Hi - l)u^^Hi)\ < di{K) + \u^^\i - l)u^^\i)\ 

(i:Re(pi(i)i)<0,pi(i)7^0) 



(i:Re(pi(i)i)>0) 



where 



With 



di{K) = max \u^'^'^ (i — l,y, e)u^'^\i,y,e) 

ie[nin + l,nfi],yeK,ee[0,eK] 



d2{K)= sup \t-t+{v,e)\'/''\u^^\t-e,v,e)u^''\t,v,e)\, 

which from the asymptotic properites of u^^\ j = 1,2 and Theorem 3.12 is finite. We find 
using the integral test that, 

d2{K)e^ 



y: \u^'\i-i)u^'^{i)\< Y 



3 , 



|t,-t+(2/,e)|V2 



(i:Re(pi (i) ^)>0) (i:Re(pi (i) ^)>0) 

< e-^d,{K) / < e-^d(K). 

A similar argument bounds the remaining sum which gives the result for j = 2. The 
argument for j = 1 is similar. □ 
We now complete the proof of Theorem A3. The bounds above, the nonvanishing, 
and asymptotic properties of u^^^ when used in equation (A. 37) imply 

n/i-l 

\(T2{n)\<d3{K)e + d4K)e'' Y ^sWkaWI. 

i=n+l 

The discrete version of Gronwall's inequality and the fact that ^[2^^-^ = c'"2(^) gives the 
bound (A. 32) for i = 2. The uniformitiy of the bound on K as well as the smoothness 
assertions follow from Theorems 3.12 and A2. The result for i = 1 follows in a similar 
manner. The proof for y e -£'2^(0) follows as above. □ 
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